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Conversion  Factors,  Non-SI  to  Sf 
Units  of  Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
units  as  follows: 


Multiply 

By 

To  Obtain 

degrees  (angle) 

0.01745329 

radians 
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1  Introduction 


Background 

Passive  classification  and  location  of  military  targets  have  significant 
advantages,  when  compared  to  active  methods  such  as  radars,  provided 
passive  techniques  can  achieve  acceptable  results.  Active  signal  tech¬ 
niques  are  vulnerable  to  interception,  “jamming,”  as  well  as  location  and 
destruction  of  the  transmitter  site.  Furthermore,  the  creation  of  an  energy 
field  requires  auxiliary  power  and  other  associated  complexities.  In  con¬ 
trast,  passive  systems  sense  only  that  energy  generated  by  the  source  itself 
(e.g.,  sound,  heat)  or  by  the  natural  environment  (e.g.,  light)  such  that 
monitoring  of  a  target  may  go  undetected.  Acoustic  techniques  have  ad¬ 
vantages  over  other  techniques,  such  as  radars,  when  terrain  or  foliage 
masking  degrades  line-of-sight  performance. 

Helicopters  represent  a  very  signiHcant  threat  to  ground  forces.  Mis¬ 
sions  can  be  conducted  at  a  low  altitude  where  radar  systems  are  ineffec¬ 
tive.  However,  helicopter  operations  are  particularly  vulnerable  to 
detection  and  potential  mine  activation,  since  they  are  often  restricted  to 
takeoff/landing  zones  and  low  altitude  avenues  of  approach. 


Purpose 

The  purposes  of  the  study  ^  presented  herein  are  as  follows: 
a.  To  present  a  method,  using  a  single  acoustic  sensor,  for 
(1)  Velocity  determination  of  approaching  helicopters. 


This  report  is  an  expansion  of  a  presentation  to  a  NATO  Conference  in  1990  (Olson 
and  Cress  1990). 
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(2)  Expected  range  at  the  closest  point  of  approach  (CPA).* 
b.  To  expand  the  characterization  of  helicopter  signatures. 


Approach 

The  approach  given  here  exploits  the  repetitive  nature  of  the  helicopter 
signature  and  uses  reasonable  assumptions  concerning  the  flight  proHle  to 
estimate  range  from  a  single  sensor.  The  technique  uses  commonly  avail¬ 
able  signal  processing  software  and  hardware,  such  as  the  Fast  Fourier 
Transform  (FFT),  to  extract  the  critical  information. 

The  approach  is  based  on  the  fact  that  any  acoustic  source  moving  in  a 
straight  line  displays  the  typical  Doppler  shift  of  sound  waves.  Figure  1 
simulates  the  frequency  of  the  main  rotor  blades  of  a  helicopter  moving  at 
33. S  m/sec  and  at  an  offset  range  of  182  m  as  it  passes  by  an  acoustic 
receiver.  As  will  be  illustrated  in  Chapter  2,  the  frequency  of  the  ap¬ 
proaching  helicopter  indicates  its  velocity.  Moreover,  the  slope  of  the 
Doppler-shifted  frequency  illustrated  in  Figure  1  near  CPA  can  be  used  to 
estimate  the  range. 

The  phenomenon  that  permits  signiEcant  improvement  in  the  calcula¬ 
tion  accuracy  is  that  of  temporal  coherence  —  some  sources  output  a  very 
coherent  (consistent  or  unchanging)  basic  frequency^  over  a  long  period 
of  time.  Evidence  indicates  that  helicopter  signals  have  long  periods  of 
temporal  coherence  (Cress  and  Olson  1990).  Changes  of  the  phase 
through  time  will  be  analyzed  in  later  sections  and  used  to  predict  the  in¬ 
stantaneous  or  short-term  frequency  to  a  high  degree  of  accuracy.  When 
applied  to  real  field  data  this  approach  can  be  used  to  determine  the  funda¬ 
mental  frequency,  velocity,  and  range  at  CPA. 


Method 

The  analysis  of  the  change  of  the  phase  through  time,  which  will  be  de¬ 
scribed  herein,  can  be  implemented  using  conventional  methods  and  hardware. 


For  convenience,  symbols  and  abbreviations  are  listed  in  the  Notation  (Appendix  D). 

2 

The  phrase  “basic  frequency”  means  the  frequency  without  any  Doppler  shift  This 
could  be  any  distinct  bud  of  energy  (e.g„  fundamental  frequency  of  main  rotor  blades, 
second  harmonic  of  main  rotor  blades,  fundamental  of  tail  rotor  blades,  etc.).  The  word 
“basic"  is  used  rather  thu  “fundamental”  to  emphasize  that  the  technique  also  worts  on 
the  higher  harmonics.  Likewise  “basic”  is  used  rather  than  “harmonics”  to  distinquish  the 
frequency  before  Doppler  shift  from  the  value  after  shifting.  (See  Appendix  A  for  techni¬ 
cal  defmitions  and  Appendix  B  for  ciariTication  of  these  frequencies.) 
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but  is  subject  to  source  parameters  (see  “Assumptions”  below).  For  this 
study,  helicopters  with  simple  flight  profiles  are  the  applicable  sources.  A 
detailed  description  of  the  equipment  (hardware)  used  in  the  field  applica¬ 
tion  is  listed  in  Table  4. 

The  theoretical  description  of  the  method  is  presented  in  Chapter  2. 

The  application  of  this  method  on  the  fundamental  of  the  main  rotor  fre¬ 
quency  of  the  moving  helicopter  is  covered  in  detail  in  Chapters  3,  4,  and 
S.  The  application  of  this  method  on  both  the  second  harmonic  of  the 
main  rotor  and  the  fundamental  of  the  tail  rotor  is  briefly  discussed  in  Ap¬ 
pendix  C. 


Issues  and  Constraints 


Assumptions 

The  velocity,  range,  and  time  of  CPA  can  be  derived  from  one  or  two 
receivers  under  certain  assumptions: 

a.  The  helicopter  is  traveling  in  a  straight  line. 

b.  The  helicopter  is  traveling  at  a  constant  velocity. 

c.  The  helicopter  possesses  a  near  constant  basic  frequency. 

One  receiver  is  sufficient  in  some  circumstances  (see  “Number  of 
receivers  required,”  page  4)  when  enough  acoustic  data  are  available  to 
calculate  all  the  information. 


Data  requirements 

If  the  basic  frequency  of  the  helicopter  is  known  from  sound  classifica¬ 
tion  (e.g.,  amplitude  versus  frequency  spectrum),  the  acoustic  information 
can  be  processed  to  determine  the  approaching  velocity  long  before  the 
source  reaches  the  receiver;  the  range  (and  time)  at  CPA  can  be  estimated 
shortly  before  the  helicopter  reaches  the  receiver.  This  information,  com¬ 
bined  with  other  acoustic  techniques  indicating  the  direction  from  which 
the  sound  is  coming  (Stabler,  Baylot,  and  Cress  1976;  Dudgeon  and  Mer- 
sereau  1984;  Smith  1984),  allows  timely  location  of  the  aircraft. 

If  the  basic  rotor  frequency  is  not  known,  the  entire  problem  can  be 
solved  by  having  two  spatially  separated  acoustic  receivers.  The  receiver 
nearest  the  approaching  aircraft  would  then  have  the  task  of  calculating 
the  basic  frequency  and  of  transmitting  this  information  back  to  the 
second  receiver.  Hence,  the  second  receiver  could  calculate  the  range  at 
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Number  of  receivers  required 


One  receiver  is  sufficient  under  either  of  the  following  assumptions: 

a.  The  basic  rotor  frequency  of  the  helicopter  is  known. 

b.  Calculation  of  the  range  and  velocity  after  the  helicopter  passes 
CPA  is  acceptable. 


Computational  accuracy 

The  formulas  to  calculate  the  velocity  and  basic  frequency  are  simple. 
Only  the  knowledge  of  the  approaching  and  retreating  Doppler-shifted  fre¬ 
quencies  and  the  velocity  of  sound  are  required.  The  range  formula  also 
requires  the  slope  of  the  frequency  at  CPA.  All  of  these  formulas  can  be 
easily  derived;  however,  they  are  very  sensitive  to  errors  in  the  Doppler- 
shifted  frequencies. 

The  primary  shortcoming  in  using  Fourier  Transform  methods  is  the  ac¬ 
curate  calculation  of  the  required  frequencies.  To  obtain  adequate  resolu¬ 
tion  of  frequency  for  the  application  presented  herein,  conventional  FFT 
techniques  would  require  observation  of  the  signature  over  an  unac¬ 
ceptably  long  time  interval.  The  time  interval  is  unacceptable  because  the 
frequency  is  constantly  changing.  Since  the  source  of  the  acoustic  signal 
in  this  investigation  is  moving,  frequencies  detected  by  a  non-moving 
receiver  are  changing  with  time  (i.e.,  a  Doppler  shift  occurs  (Weidner  and 
Sells  1965)).  Such  a  shift  in  the  frequencies  causes  a  smearing  by  the 
FFT  of  the  power  into  a  number  of  different  frequency  bins,*  decreasing 
the  accuracy  for  conventional  calculations.  Such  smearing  is  unaccep¬ 
table  for  this  implementation. 


Recording  accuracy 

The  addition  of  some  noise  onto  the  recorded  signal  impairs  an  ac¬ 
curate  calculation  of  the  frequency  at  any  time  using  the  standard  FFT 
techniques.  Wind  can  be  a  major  source  of  noise  and  Doppler  smearing 
by  moving  media.  The  presence  of  wind  may  introduce  two  impediments: 
a  deterioration  of  the  signal-to-noise  ratio  and  changes  in  the  travel  time 
and  frequency  (Doppler  shift)  of  the  sound.  These  complications  are  not 
considered  in  this  report.  A  poor  signal-to-noise  ratio  is  not  expected  to 
be  a  serious  problem  for  ranges  less  than  500  m.  However,  the  wind  effect 


The  concept  of  “frequency  bins”  is  the  discrete  frequency  classification  that  the  FFT 
assigns  to  a  real  signal.  For  example,  a  1 -sec-long,  19.1 -Hz  signal  would  have  most  of  its 
energy  assigned  to  the  19-Hz  FFT  bin  with  the  second-largest  amount  of  energy  assigned 
to  the  20-Hz  FFT  bin. 
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on  acoustic  velocity  could  be  a  problem,  particularly  when  the  wind  speed 
is  greater  than  10  percent  of  the  aircraft  velocity. 


Potential  Applications 

Potential  applications  include  reconnaissance  sensors  assigned  to  per¬ 
form  early  detection  and  analysis  of  a  helicopter  acoustic  signal.  Another 
use  would  be  for  sensors  assigned  the  decision  of  whether  or  not  to  ac¬ 
tivate  mines.  Potential  locations  include  situating  the  sensors  in  a  valley 
ti.ai  serves  as  a  natural  avenue  of  approach  below  radar  or  near  likely  land¬ 
ing  zones. 
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2  Analytical  Method 


The  mathematical  foundation  for  the  detailed  analysis  includes  an  ele¬ 
mentary  review  of  the  Doppler  shift  theory  of  acoustic  waves.  The  for¬ 
mulas  (and  also  the  time  when  enough  acoustic  data  has  been  collected  to 
perform  the  necessary  calculations)  depend  on  whether  or  not  the  type  of 
helicopter  (and  its  attendant  acoustic  signal)  is  known.  Both  cases  require 
an  accurate  knowledge  of  the  Doppler-shifted  frequencies. 

To  obtain  the  high  degree  of  accuracy  in  the  frequency  estimate,  this  ap¬ 
proach  exploits  conventional  signal-processing  theory.  A  close  investiga¬ 
tion  of  the  phase  information  derived  by  the  FFT  of  a  constant-frequency 
wave  provides  the  basis  for  the  significantly  improved  accuracy. 


Doppler  Theory 

The  most  general  case,  the  non-recognized  helicopter,  requires  very  lit¬ 
tle  knowledge  of  the  approaching  target.  Two  spatially  separated 
receivers  are  suited  for  this  situation.  Other  alternatives  are  use  of  a 
single  receiver  with  correct  classification  of  the  helicopter  using  the  acous¬ 
tic  signature  or  use  of  a  supplemental  sensor,  such  as  a  thermal  infrared 
(IR)  radiometer  and  determination  of  the  acoustic  fundamental  from  the 
modulation  of  the  blade  rotation.  If  two  receivers  are  used,  the  front 
receiver  has  the  task  of  calculating  the  basic  frequency  of  the  helicopter 
and  relaying  it  to  the  second  acoustic  receiver  for  additional  calculations. 
Both  receivers  may  also  be  able  to  independently  calculate  the  velocity, 
range  at  CPA,  and  time  of  CPA.  For  the  approach  presented  herein,  the 
front  receiver  will  not  be  able  to  calculate  Ais  information  until  after  the 
helicopter  has  passed.  The  second  receiver,  in  the  unknown  helicopter 
case,  can  calculate  this  information  just  before  the  helicopter  reaches 
CPA.  This  second  receiver  can  then  be  programmed  for  the  case  where 
the  type  of  helicopter  is  known  and  only  one  microphone  is  required. 
Hence,  only  one  receiver  will  be  considered  here. 
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Assumptions 


Table  1  summarizes  the  assumptions  required  for  the  following  calcu¬ 
lations.  The  principal  reason  for  assuming  only  one  moving  source  is  to 
ensure  the  possibility  of  separating  the  signals  in  the  frequency  domain. 
Extension  to  multiple  sources  at  the  same  acoustic  generation  frequencies 
will  require  further  study.  This  separation  is  needed  to  accurately  deter¬ 
mine  the  approaching  and  retreating  far-field  frequencies  and  the  slope  of 
the  frequency  near  CPA.  The  second  assumption  of  a  constant  velocity 
and  direction  is  probably  not  valid  over  long  periods  of  time  but  is  ap¬ 
plicable  for  short  intervals.  An  analysis  of  the  data  can  usually  determine 
whether  this  assumption  is  valid. 


• 

Table  1 

Assumptions  for  One  Acoustic  Receiver 

B 

There  is  only  one  moving  noise  source. 

B 

The  helicopter  is  traveling  at  a  constant  velocity  (V)  in  a  straight  line. 

B 

The  velocity  of  sound  (C)  from  the  helicopter  to  the  receiver  is  constant  and  is  known. 

B 

V  is  less  than  C. 

B 

The  helicopter  outputs  a  basic  frequency  (fo)  that  remains  constant  when  measured  by 
an  observer  in  the  helicopter  (i.e.,  the  helicopter  point  of  view). 

The  last  assumption  in  Table  1  implies  that  any  frequency  change  mea¬ 
sured  at  a  fixed  receiver  is  caused  by  a  Doppler  shift.  The  validity  of  this 
assumption  for  a  Chinook  helicopter  (CH-46)  was  one  of  the  main  conclu¬ 
sions  of  Cress  and  Olson  (1990):  that  the  source  emits  a  temporally  coher¬ 
ent  signal  at  the  fundamental  over  time  periods  on  the  order  of  IS  sec. 


Calculations 

Figure  2  is  derived  from  a  model  subject  to  the  assumptions  in  Table  1. 
In  this  application,  the  term  “far  field”  means  that  the  separation  between 
the  source  and  the  receiver  is  enough  to  ensure  that  the  Doppler-shifted 
frequency  is  less  than  10  percent  of  the  total  change  of  the  Doppler  shift 
for  the  approaching  frequencies  relative  to  the  retreating  frequencies.  In 
both  cases  the  receiver  records  a  signal  with  essentially  a  constant  fre¬ 
quency.  Standard  Doppler  theory  yields  the  following  formulas: 


(1) 


(2) 
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where 


fa  =  the  approaching  frequency 
fa  -  basic  acoustic  frequency  of  the  source 
V  =  velocity  of  that  source 
C  =  velocity  of  sound  in  still  air 
fr  =  retreating  frequency 

These  two  formulas  can  be  solved  for  and  V  in  terms  of/^  ,/^ ,  and  C: 


^  _1AA 

fa^fr 


V  = 


fa-fr 

fa^fr 


(3) 

(4) 


The  aircraft  at  CPA  is  sending  out  a  signal  to  the  receiver  at  frequency 
fa.  A  more  involved  calculation  subjected  to  the  assumption  in  Table  1 
shows  that  when  the  helicopter  passes  the  receiver  at  a  range  R  at  CPA, 
the  slope  of  the  change  in  the  frequency  /at  CPA  determines  the  range  (R): 


R  = 


c  •  f  (CPA) 


(5) 


Time  delay  from  aircraft  to  receiver 

The  time  required  for  sound  to  travel  from  the  helicopter  to  the 
receiver  actually  introduces  a  slight  change  in  the  shift  of  the  frequencies. 
However,  this  shift  has  no  impact  on  the  far-field  approaching  and  retreat¬ 
ing  frequencies,  and  minimal  impact  on  the  slope  of  the  Doppler  shift  at 
CPA.  This  difference  is  illustrated  in  Figure  3  for  the  situation  shown  in 
the  acoustic  plot  of  Figure  1.  Consequently,  the  time  shift  will  be  ignored 
in  the  derivation. 


Calculations  when  the  basic  frequency  Is  known 

As  soon  as  the  far-field  approaching  frequency  fg  is  detected  and  ana¬ 
lyzed,  the  approaching  velocity  V  can  be  calculated  with  the  following 
formula: 
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This  formula  is  derived  from  Equation  1.  The  formula  for  the  range  at 
CPA  remains  unchanged  from  Equation  5  and  is  repeated  below: 

-  f  . 

_  JO _ 1 - 

C  •  ^  {CPA) 


Time  of  CPA 


It  should  be  noted  that  when  the  aircraft  is  at  CPA,  it  is  sending  a  fre¬ 
quency  to  the  receiver  that  is  not  frequency  shifted.  If  represents  the 
time  at  which  the  ground  is  receiving  a  Doppler-shifted  frequency 
and  if  equals  then  the  sound  was  generated  when  the  aircraft  was 
at  CPA.  Since  that  sound  travels  a  distance  /?  at  a  velocity  C,  the  actual 
time  the  aircraft  was  at  CPA  (t^PA) 


*CPA 


\  J 


(7) 


If  the  slope  of  the  frequency  remains  approximately  constant  for  some 
time  before  CPA,  the  actual  time  and  range  at  CPA  of  the  aircraft  can  be 
estimated  before  it  reaches  that  point  using  the  slope  from  the  earlier 
times  (t): 


m)  = 


(8) 


The  slope  of  the  line  tangent  to  the  current  frequency  curve  at  a  time  t 
before  CPA  can  also  be  used  to  estimate  the  slope  between  the  current 
point  {t,f(t))  on  the  curve  and  the  point  {tg,!^  of  CPA  (see  Figure  4): 


4L 

dt 


(0 


fo-m 

to-f 


(9) 


Solving  Equation  9  for  yields  the  estimate  tj{t)  (made  at  time  t)  of  when 
the  sound  made  at  CPA  will  reach  the  receiver 


fo-m 
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Using  this  value  and  Equation  8  in  the  formula  for  the  actual  time  that 
the  aircraft  is  at  CPA,  Equation  7,  yields  the  estimate  at  time  t  of  the  time 


of  CPA 


=  t  + 


fo-m 


-fo  • 

^  ■  f « 


which  can  be  simplified  to 


^CPA  (0  =  ^  + 


\fo -m]  •  +fo  ■ 

C’  •  f  (') 


(11) 


This  equation  is  the  one  of  interest  if  the  intention  is  to  intercept  the 
aircraft.  It  represents  the  estimate  (made  at  a  time  t)  of  when  the  aircraft 
will  actually  arrive  at  CPA.  The  key  assumption  of  this  formula  is  that  the 
slope  of  the  frequency  curve  approximates  the  corresponding  slope  when 
the  aircraft  is  at  CPA.  The  expected  range  is  given  in  Equation  8. 


Accurate  Calculation  of  Frequency 

The  previous  formulas  are  of  limited  value  if  the  frequency  is  not 
known  accurately.  The  standard  FFT  calculation  over  an  interval  of  time 
has  a  relationship  that  limits  the  resolution  of  the  frequency: 

1  =  [Frequency  Increment  (Hz)]  •  [Time  Interval  Length  {sec)]  (12) 

This  formula  suggests  that  resolving  a  signal  to  within  0.1  Hz  requires 
applying  the  FFT  to  a  10-sec-long  signal.  This  may  be  possible  in  the  far- 
fleld  situation,  but  it  is  definitely  not  applicable  near  CPA  for  a  moving 
helicopter  on  the  order  of  30  m/sec  or  greater  at  ranges  of  several  hundred 
meters  or  less.  Near  CPA  the  required  accuracy  on  the  frequency  is  much 
higher  since  the  slope  of  the  frequency  must  be  precisely  calculated  to 
determine  the  range  (Equation  8)  and  time  of  CPA  (Equation  11).  Unfor¬ 
tunately,  near  CPA  the  frequency  of  the  signal  is  changing  rapidly;  thus, 
the  FFT  of  a  long  time  interval  will  represent  only  an  average  of  that  time 
interval  and  will  deEnitely  not  represent  an  estimate  of  the  instantaneous 
frequency  required  for  the  above  formulas. 

Selection  of  a  short  time  interval  for  the  FFT  results  in  poor  resolution 
of  the  frequency  if  the  conventional  calculations  are  assumed  (see  Figure  5). 
If  the  time  interval  is  only  O.S  sec  in  length,  the  resolution  of  the  frequen¬ 
cy  provided  by  the  FFT  is  only  2  Hz.  It  should  be  noted  that  the  entire 
Ooppler  shift  of  the  strongest  signal  (Figure  1)  is  less  than  5  Hz.  Clearly, 
a  2-Hz  resolution  of  the  frequency  is  useless  for  calculating  the  slope  at 
CPA  for  that  signal. 
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An  improved  frequency  resolution  can  be  obtained  with  the  concepts  of 
the  temporal  coherence  of  the  source  signal  and  a  careful  analysis  of  the 
phase  of  the  cross-spectral  density  (CSD).  A  simple  example  illustrates 
these  concepts.  The  method  exploits  the  assumption  that  the  bandwidth  of 
the  unknown  source  basic  frequency  is  significantly  narrower  than  the  fre¬ 
quency  lines  in  the  FFT. 


Examination  of  the  FFT 


FFT  of  a  19.1-Hz  signal.  Suppose  the  FFT  was  used  to  analyze  a  sig¬ 
nal  h(t)  which  is  actually  a  pure  19.1 -Hz  cosine  wave: 

h{f)  =  cos  (271  •  19.1  •  t) 


Assume  that  successive  1-sec  intervals  will  be  used  to  analyze  this  signal. 
According  to  Equation  12,  this  interval  length  permits  a  1-Hz  frequency 
resolution.  Consider  the  19-Hz  term  of  the  FFT.  The  FFT  will  identify 
two  real  numbers  A  and  6  for  which  the  function  g(t)  closely  matches  h(t). 

g(t)  =  A  cos  (271  •  19  •  /  +  0) 


Significance  of  phase  term.  Clearly  the  amplitude  A  should  be  ap¬ 
proximately  1,  but  0  will  definitely  depend  on  the  time  interval  [n,n+l] 
over  which  the  FFT  was  calculated  (see  Figure  6).  By  rewriting  h(t)  as 


h{t)  =  cos 


271  •  19  •  7  +  271 


j_ 

10 


(13) 


V  / 

it  seems  reasonable  that  0  represents  the  average  of  the  term  27i-t/10  over 
the  time  interval  [n,n+l] 


Bin)  = 


271 


10 


27t  •  (w  +  1) 


10 


(14) 


=  7t 


(2n  -H  1) 
10 


Graphically  this  implies  that  the  FFT  frequency  domain  transformation 
of  successive  time  intervals  rotates  around  the  origin  of  the  complex  plane 
at  a  constant  rate.  The  first  point  plots  at  an  angle  of  7i/10  rad  with 
respect  to  the  positive  x-axis.  The  succeeding  points  plot  at  37i/10,  Sti/IO, 
77t/10.  971/10,  ll7t/10,  1371/10,  1571/10.  177t/10,  1971/10,  Ti/lO,  37t/10, 

571/10, ...  for  the  [n,n+l],  n=  1,  2, . . .  time  intervals.  It  takes  10  sec  for 
the  19-Hz  FFT  frequency  domain  value  to  move  around  the  complex  unit 
circle  (Ahlfors  1966)  when  analyzing  the  19.1-Hz  signal  (see  Figure  7). 
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Impact  on  the  CSO 


Phase  changes  with  time.  This  constant  phase  rotation  of  the  angles 
of  the  FFT  response  over  succeeding  time  intervals  indicates  that  the  dif¬ 
ference  ^  in  these  phases  remains  constant.  For  example,  for  the  19.1 -Hz 
signal,  the  difference  between  the  total  phase  estimates  for  succeeding 
time  intervals  n  to  n+l  is  given  by: 


(15) 


Now  Equations  13  and  14  can  be  written  as 

0(/i)  =  0(0)  +  n  •  p  (16) 

and 

h{t)  =  cos  (27c  •  19  •  f  +  P  •  t) 


This  formula  appears  to  be  close  to  the  original  definition  of  h(t)  in 
Equation  13.  However,  Equation  17,  which  expresses  a  completely  differ¬ 
ent  point  of  view,  states  that  if  the  phase  differences  of  the  19-Hz  FFT 
of  successive  time  intervals  are  constantly  equal  to  p  radians,  then  a  (19 
+  p/2n)-Hz  signal  might  exist  in  the  time  domain. 

Definition  of  CSD.  If  //^(f)  is  deflned  to  be  the  complex  number  that 
is  the  output  of  the  FFT  applied  to  the  signal  h(t)  at  f  Hz  of  the  time  inter¬ 
val  [nji+l],  then  the  CSD  of  this  signal  and  frequency  over  that  time  inter¬ 
val,  is  deHned  as 

^rff)  =  Hnif)  •  d*) 

where  the  *  denotes  the  complex  conjugate.  The  phase  of  the  CSD  repre¬ 
sents  the  difference  in  the  phases  of  successive  time  intervals,  [n-l,n]  to 
[njt+l],  of  the  FFT  of  the  frequency  f.  This  suggests  that  the  phase  of 
represents  the  instantaneous  change  of  the  phase  of  the  f-Hz  com¬ 
ponent  of  the  h(t)  signal  at  time  t  =  n. 

In  the  example  of  the  signal  h(t)  representing  the  19.1-Hz  cosine  wave, 
the  CSD  of  19  Hz,  Xn(19)  has  a  constant  phase  of 


P  = 


2k 

10 


radians 


=  0.1  cycles 
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Likewise,  if  a  pure  (19+5f)-Hz  signal  with  I  5f  l<.5  is  transformed  to 
the  frequency  domain  and  analyzed  with  the  CSD  of  adjacent  1-sec  inter¬ 
vals,  then  the  phase  of  the  19-Hz  CSD  will  be  constant  at  5f  cycles. 

Complex  plane  graph  over  time  of  the  CSD.  Since  the  plot  of  the 
FFT  rotates  around  the  unit  circle  at  constant  phase  increments  and  the 
CSD  represents  the  phase  difference  of  the  FI^,  logically  the  CSD  will  be 
a  constant  for  the  pure  signal.  Hence,  all  the  CSD  points,  X„(19)  n=l,  2, . . . 
of  a  pure  (19-i-5f)-Hz  signal  in  the  complex  plane  will  actually  plot  on  top 
of  each  other  5/  of  the  way  around  the  unit  circle  (see  Figure  7)  if  the 
amplitude  does  not  change  with  time.  In  the  presence  of  a  small  amount 
of  recording  noise  or  possibly  a  slightly  impure  signal,  all  these  points 
will  plot  in  roughly  the  same  location  of  the  complex  plane. 

Cumulative  sum  of  the  CSD.  In  the  presence  of  moderate  noise,  the 
CSD  signal  would  plot  over  a  much  larger  range  of  angles  in  the  complex 
plane.  One  way  to  determine  the  overall  trend  is  to  use  complex  addition 
over  N+l  adjacent  time  intervals  to  add  the  signal  and  cancel  the  noise: 

N  N  (19) 

Rs)  n=\ 


The  plot  of  the  cumulative  sum  of  the  CSD  (Figure  8)  will  clearly  show 
a  trend  for  a  narrow  band  frequency  input  with  a  good  signal-to-noise 
ratio.  Successive  points  should  plot  in  roughly  a  straight  line  radiating 
away  from  the  origin  (Figure  8).  That  figure  plots  the  locations  of  the 
sum  for  N=l,  2, . . .  14. 


Important  inverse  problem.  Up  to  this  point  the  frequency  content  of 
the  time  signal  was  usually  known.  Suppose  the  current  task  was  to  ana¬ 
lyze  the  frequency  content  of  an  unknown  signal  h(t).  If  the  FFT  is  calcu¬ 
lated  for  two  successive  time  intervals  [n-l,n]  and  [n,n-i-l]  and  the  largest 
amplitude  stays  at  19  Hz,  can  a  better  estimate  of  the  frequency  be  made? 

The  above  calculations  indicate  that,  when  measured  in  cycles,  the 
phase  of  the  CSD  of  the  19-Hz  frequency  between  these  two  adjacent  time 
intervals  is  really  the  difference  between  the  actual  frequency  content  of 
the  signal  (measured  in  Hz)  and  19  Hz: 

If  19  Hz  is  the  largest  CSD  amplitude  of  a  real  signal,  then 


Actual  Frequency  {Hz)  =  \9Hz  + 


CSD  Phase  of  19  Hz  {cycled 
Length  of  FFT  Interval  (sec) 


(20) 


where  the  phase  of  the  CSD  is  between  -.5  and  .5. 


For  a  signal  with  good  signal-to-noise  ratio,  the  above  computation  can 
be  made  for  each  successive  CSD  value  (between  each  successive  FFT 
interval).  The  presence  of  noise  does  degrade  the  estimate  of  the  actual 
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frequency.  Assuming  that  the  signal  is  not  changing  significantly  over 
time,  the  cumulative  sum  of  a  number  of  CSD  values  should  have  a  phase 
that  minimizes  the  corruption  of  the  noise,  yielding  a  better  phase  value 
than  the  individual  CSD  numbers  in  Equation  20. 


Temporal  coherence  implications 


The  sum  of  the  products  of  the  FFT  of  a  signal  h(t)  multiplied  by  the 
complex  conjugate  of  the  preceding  interval  (20)  actually  represents  the 
numerator  of  the  temporal  coherence,  TemCoh(0.  of  the  signal  h(t)  at  the 
frequency  f 


TemCoh{f) 


N 

X 

ffsl _ _  _ 


Vi 

'  N 

• 

/r=l 

n=l 

(21) 


The  previous  calculations  indicate  that  if  a  signal  has  a  consistent  fre¬ 
quency  over  time  close  to  f,  then  the  temporal  coherence  should  have  a 
constant  phase  and  an  amplitude  close  to  one.  In  fact,  that  phase  mea¬ 
sured  in  cycles  defines  how  much  the  original  signal’s  frequency  (mea¬ 
sured  in  hertz)  differs  from  f. 

On  the  other  hand,  if  the  amplitude  of  the  temporal  coherence  is  not 
close  to  unity,  then  neither  the  phase  of  the  temporal  coherence  nor  the 
phase  of  the  cumulative  sum  of  the  CSD  over  that  long  a  collection  of 
time  intervals  should  be  used  as  a  replacement  for  the  CSD  phase  in  Equa¬ 
tion  20.  A  plot  of  the  graph  of  the  cumulative  sum  of  the  CSD  (e.g..  Fig¬ 
ure  9)  shows  good  coherence  when  the  plot  radiates  in  approximately  a 
straight  line  from  the  origin. 


Amplitude  significance 

It  is  important  to  include  the  amplitude  in  the  analysis  along  with  phase 
and  temporal  coherence.  The  significance  can  be  directly  seen  by  con¬ 
sidering  the  pure  19.1 -Hz  cosine  signal  again. 

The  20-Hz  FFT  analysis  of  the  19.1-Hz  signal.  The  same  analysis 
can  be  performed  on  the  same  signal — but  using  the  20-Hz  FFT  and  CSD 
analysis.  Our  time  signal  h(t)  can  be  rewritten  so  that  the  20  Hz  will  see 
it: 
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hit)  =  cos(2n  ■  19.1  ■  t) 


(22) 


=  cos  (2ji  •  20  •  f  -  2ji  ■  0.9  ■  t) 

Hence,  the  FFT  will  see  a  constant  phase  shift  of  -271-0.9  rad  (or 
-0.9  cycles)  over  successive  1 -sec-long  intervals.  Again  the  phase  change 
of  the  FFT  is  the  phase  of  the  CSD.  Since  the  CSD  and  FFT  phases  are 
only  analyzed  to  within  one  revolution  of  the  unit  circle,  this  phase 
change  of  the  FFT  will  be  reported  by  the  CSD  to  be  +27t  •  0.1  rad  (or 
+0.1  cycle).  This  value  for  20  Hz  is  exactly  the  same  CSD  phase  seen  in 
the  19-Hz  CSD  analysis. 

Since  the  CSD  for  20  Hz  is  constant  at  0.1  cycle  for  all  1 -sec -long  time 
intervals,  the  cumulative  sum  of  all  CSD  values  will  align,  and  the  tem¬ 
poral  coherence  will  be  1.0  also.  One  might  be  tempted  to  claim  (in  the 
spirit  of  Equation  20  if  the  amplitude  was  ignored)  that  we  are  actually 
analyzing  a  pure  20.1 -Hz  signal.  The  major  difference  between  the  FFT, 
CSD,  and  cumulative  sum  of  the  CSD  for  20  Hz,  and  their  counterparts 
for  19  Hz,  is  the  amplitude  of  all  these  complex  numbers.  The  19-Hz 
amplitudes  will  be  signiffcantly  larger. 

CSD  phase  for  all  frequencies.  In  theory,  the  CSD  phase  of  our 
19.1-Hz  signal  will  be  exactly  the  same  for  all  1-Hz  frequency  analysis 
values  in  the  frequency  domain  (e.g.,  the  18-Hz  and  21 -Hz  CSD  phase 
will  also  be  0. 1  cycle).  Likewise  the  phase  of  the  cumulative  sum  will 
always  be  0.1. 

The  words  “in  theory”  preface  this  section  because,  in  reality,  the 
machine  significance  restrictions  of  all  computers  will  show  some  dif¬ 
ferences  in  the  phase.  The  difference  should  increase  as  the  analysis  fre¬ 
quency  moves  away  from  the  pure  19.1-Hz  signal  (see  Table  2). 

CSD  amplitude  for  all  frequencies.  The  only  difference  in  these  com¬ 
plex  numbers  is  in  their  amplitude.  The  closer  the  analysis  frequency  is  to 
19.1  Hz,  the  higher  the  amplitude  of  the  complex  number  will  be.  See 
Table  3  for  a  summary  of  the  coherence  amplitude  (over  fifteen  1 -sec-long 
intervals)  and  the  first  few  CSD  amplitudes  and  phases  for  the  pure 
19.1-Hz  cosine  signal. 

Amplitude  of  TemCob  for  all  frequencies.  Finally,  the  amplitude  of 
the  temporal  coherence  is  useless  (in  the  case  of  a  pure  unchanging  signal) 
for  determining  which  analysis  frequency  is  important.  The  amplitude  of  the 
coherence  (the  coherence  column  of  Tables  2  and  3)  is  close  to  1.0  for  all  fre¬ 
quencies  near  19.1  Hz. 
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Table  2 

Cross-Spectral  Density  Phase  of  a  19.1-Hz  Cosine  Function 


PhaM  (dagrMs)  for  Successive  Time  Intervals  (sec) 


Frequency 

(Hx) 

Coherence 
(15  sec) 

10,1] 

11.2] 

[1.2] 

[2.3] 

[2.3] 

[3.4] 

[3.4] 

[4.5] 

[4.5] 

[5.6] 

[5.6] 

[6.7] 

10.0 

0.881 

25.4 

52.7 

55.7 

26.9 

19.3 

25.4 

11.0 

0.908 

27.1 

50.4 

52.8 

28.5 

21.2 

27.1 

12.0 

0.932 

28.7 

48.2 

50.1 

30.0 

23.1 

28.7 

13.0 

0.952 

30.1 

46.1 

47.6 

31.2 

24.9 

30.1 

14.0 

0.968 

31.4 

44.2 

45.3 

32.3 

26.8 

31.4 

15.0 

0.980 

32.5 

42.4 

43.1 

33.3 

28.6 

32.5 

16.0 

0.989 

33.5 

40.7 

41.2 

34.1 

30.5 

33.5 

17.0 

0.995 

34.4 

39.1 

39.4 

34.9 

32.3 

34.4 

18.0 

0.999 

35.2 

37.6 

37.7 

35.5 

34.1 

35.2 

19.0 

1.000 

35.9 

36.1 

36.1 

36.0 

35.8 

35.9 

20.0 

0.999 

36.5 

34.8 

34.7 

36.4 

37.6 

36.5 

21.0 

0.997 

37.1 

33.6 

33.4 

36.7 

39.3 

37.1 

22.0 

0.993 

37.5 

32.4 

32.1 

37.0 

41.0 

37.5 

23.0 

0.989 

37.9 

31.3 

30.9 

37.2 

42.7 

37.9 

24.0 

0.983 

36.2 

30.2 

29.8 

37.3 

10111 

38.2 

25.0 

0.977 

38.5 

29.3 

28.8 

37.4 

46.1 

38.5 

26.0 

0.970 

38.7 

28.3 

27.9 

37.4 

47.7 

38.7 

27.0 

0.963 

38.9 

27.5 

27.0 

37.4 

49.3 

38.9 

28.0 

0.956 

39.0 

26.6 

26.1 

37.4 

50.9 

39.0 

29.0 

0.949 

39.1 

25.8 

25.3 

37.3 

52.5 

39.1 

30.0 

0.941 

39.1 

25.1 

24.6 

37.2 

54.0 

39.1 
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Table  3 

Cross-Spectral  Density  Amplitude  of  a  19.1 -Hz  Cosine  Function 

Frequency 

(Hz) 

Coherence 
(15  see) 

Time  intervals  (sec)  Being  Multiplied 

10.1] 

11.2] 

HhM 

[2,3] 

[2.3] 

[3.4] 

[3.4] 

[4.5] 

[4.5] 

[5.6] 

[5.6] 

[6.7] 

10.0 

0.881 

0.00014 

0.00008 

0.00007 

0.00014 

0.00019 

0.00014 

11.0 

0.908 

1 

0.00018 

0.00010 

0.00010 

0.00017 

0.00022 

0.00018 

12.0 

0.932 

0.00022 

0.00014 

0.00014 

0.00021 

0.00027 

0.00022 

13.0 

0.952 

0.00029 

0.00020 

0.00020 

0.00028 

0.00035 

0.00029 

14.0 

0.968 

0.00041 

0.00031 

0.00030 

0.00040 

0.00047 

0.00041 

15.0 

0.980 

0.00062 

0.00049 

0.00049 

0.00061 

0.00070 

0.00062 

16.0 

0.989 

0.00106 

0.00090 

0.00089 

0.00105 

0.00116 

0.00106 

17.0 

0.995 

0.00227 

0.00204 

0.00203 

0.00225 

0.00241 

0.00227 

18.0 

0.999 

0.00814 

0.00770 

0.00768 

0.00810 

0.00839 

0.00814 

19.0 

1.000 

0.96903 

0.96436 

0.96407 

0.96857 

0.97164 

0.96903 

20.0 

0.999 

0.01179 

0.01229 

0.01232 

0.01183 

0.01151 

0.01179 

21.0 

0.997 

0.00261 

0.00284 

0.00286 

0.00263 

0.00248 

0.00261 

22.0 

0.993 

0.00110 

0.00126 

0.00127 

0.00112 

0.00102 

0.00110 

23.0 

0.989 

0.00060 

0.00071 

0.00072 

0.00061 

0.00055 

0.00060 

24.0 

0.983 

0.00038 

0.00046 

0.00047 

0.00039 

0.00033 

0.00038 

25.0 

0.977 

0.00026 

0.00033 

0.00033 

0.00026 

0.00022 

0.00026 

26.0 

0.970 

0.00019 

0.00025 

0.00025 

0.00019 

0.00016 

0.00019 

27.0 

0.963 

0.00014 

0.00019 

0.00020 

0.00015 

0.00012 

0.00014 

28.0 

0.956 

0.00011 

0.00015 

0.00016 

0.0001 1 

0.00009 

0.00011 

29.0 

0.949 

0.00009 

0.00013 

0.00013 

0.00009 

0.00007 

0.00009 

30.0 

0.941 

0.00007 

0.00011 

0.00011 

0.00008 

0.00006 

0.00007 
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Analytical  Summary 


Using  the  temporal  coherence  and  a  detailed  analysis  of  the  signal 
phase  permits  a  more  detailed  resolution  of  the  frequency  of  a  signal.  The 
approach  presented  here  and  for  the  subsequent  measured  data  is  to  per¬ 
form  an  FFT  over  1 -sec-long  time  intervals.  Next,  the  cross-spectral  den¬ 
sity  of  adjacent  intervals  is  calculated.  The  signal  of  interest  is  located  by 
choosing  that  integer  frequency  yielding  the  peak  amplitude  of  the  cross- 
spectral  density.  Then  the  phase  of  the  CSD  of  that  signal  is  used  to  ob¬ 
tain  a  resolution  of  the  frequency  beyond  the  1-Hz  limitation  of  the  FFT. 

If  the  signal  remains  coherent  over  a  length  of  time  (as  determined  by  the 
temporal  coherence),  the  amplitude  and  phase  of  the  cumulative  sum  of 
the  CSD  yields  an  improved  estimate. 

In  theory,  the  resolution  is  infinite  if  the  source  signature  is  infmitely 
narrow.  In  practice,  the  noise  on  the  data,  the  resolution  of  the  digitizing 
process,  the  bandwidth  of  the  source,  and  the  computer  machine  signifi¬ 
cance  will  limit  the  frequency  accuracy.  The  application  of  this  approach 
to  the  measured  acoustic  data  that  follow  shows  improvement  over  conven¬ 
tional  processing  even  with  the  corruption  of  the  data  by  wind  and  a  60-Hz 
noise  source,  and  with  only  four  digits  of  accuracy  coming  from  the 
digitizing  process. 

This  improved  resolution  of  a  Doppler-shifting  frequency  permits  the 
calculation  of  the  basic  frequency,  velocity,  time  of  CPA,  and  range  at 
CPA  of  a  source  moving  in  a  straight  line  (and  therefore  the  range  as  a 
function  of  time).  All  of  these  data  can  be  calculated  with  one  receiver. 

The  basic  source  frequency  (i.e.,  not  Doppler-shifted  frequency)  and 
velocity  can  be  calculated  from  the  knowledge  of  the  Doppler-shifted  fre¬ 
quencies.  The  range  at  CPA  and  time  of  CPA  can  be  calculated  from  the 
knowledge  of  the  basic  frequency,  velocity,  and  slope  of  the  frequency 
near  CPA. 

If  the  basic  frequency  of  the  source  is  known  a  moderate  time  before  it 
reaches  CPA,  the  velocity  can  be  calculated  when  the  far-field  approach¬ 
ing  Doppler-shifted  frequency  is  known.  As  the  source  approaches  CPA, 
calculating  estimated  time  of  CPA  and  estimated  range  at  CPA  is  possible. 
These  estimates  improve  as  the  source  approaches  CPA. 
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Analysis  of  Far-Field 
Frequency  for  Measured 
Data 


Statistics  on  Vicksburg  Heiicopter  Test 

Acoustic  measurements  were  obtained  in  June  1990  by  personnel  of  the 
Battlefield  Environment  Group,  Environmental  Laboratory,  U.S.  Army  En¬ 
gineer  Waterways  Experimental  Station.  The  acoustic  data  were  recorded 
at  the  Vicksburg  Municipal  Airport  on  U.S.  Highway  61  south  of  the  city 
of  Vicksburg. 


Test  description 

Figure  10  illustrates  the  field  geometry  of  the  run.  All  of  the  infor¬ 
mation  presented  below  was  gathered  from  one  test.  The  statistics  of 
Run  21  are  presented  in  Table  4  and  include  data  on  the  recording  equip¬ 
ment,  source,  flight  path,  and  weather. 

As  previously  stated,  one  sensor  does  not  permit  determination  of  the 
direction  of  the  sound.  However,  using  the  ITT  enables  determination  of 
how  the  frequency  content  changes  with  time.  The  analysis  of  the  Dop¬ 
pler  shift  of  the  frequency  collected  from  one  sensor  permits  an  estimation 
of  the  velocity  and  range  of  the  source. 
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Table  4 

Statistics  Concerning  the  Field  Data 

Kam 

Statistic 

Location: 

Vicksburg  Municipal  Airport 

Date: 

June  21, 1990 

Time: 

1623  Central  Daylight  Time 

Test  No. 

21 

Audio  phone: 

Bruel  &  Kjaer,  Model  4921 

Channel: 

1 

Amplifier  gain: 

2 

Digitizing  rate: 

2,048  samples/sec 

Recording  length: 

135  sac 

Source: 

Helicopter  with  basic  frequency  172  Hz 

Altitude: 

152.4  m 

Horizontal  offset  at  CPA: 

100.0  m  east 

Source  velocity: 

33.5  m/sec 

Travel  direction  and  range: 

2  km  south  to  2  km  north  of  CPA 

Temperature: 

33.85  'C 

Relative  humidity: 

58.58%  (3  m  above  ground) 

Barometric  pressure: 

1,071  millibars 

Wind  velocity: 

0.641  m/sec 

Sound  velocity  (calculated): 

352.6  m/sec 

Signal  characteristics 

The  data  were  digitized  at  2,048  samples  per  second,  yielding  a  Ny> 
quist  frequency^  of  1,024  Hz.  Successive  1 -sec -long  sections  of  the  data 
were  analyzed  using  the  FFT  to  compute  the  frequency  content  for  that 
time  interval.  This  permits  a  frequency  resolution  of  1  Hz. 

A  detailed  plot  of  the  amplitude  as  a  function  of  time  interval  and  fre¬ 
quency  is  illustrated  in  Figure  11.  The  time  interval  from  the  beginning  of 
the  test  at  time  1  until  13S  sec  later  is  plotted  against  the  frequency  con¬ 
tent  from  1  to  59  Hz.  Each  grid  point  on  the  3-D  surface  represents  the 
FFT  amplitude  (using  a  linear  scJe)  for  a  1-Hz  increment  of  a  1 -sec-long 


A  frequency  associated  with  the  rate  of  digitization.  Frequencies  larger  than  the  Ny- 
quist  frequency  must  be  removed  before  digitization  (Sheriff  1984). 
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interval.  The  higher  the  point,  the  larger  the  amplitude.  The  maximum 
value  was  0.900.*  All  amplitudes  below  a  threshold  value  of  0.033  were 
set  to  zero. 

A  number  of  bands  of  energy  show  the  expected  signature  of  Doppler- 
shifted  acoustic  energy  from  a  source  traveling  in  a  straight  line.  From 
time  10  sec  to  65  sec,  the  strongest  frequency  remains  approximately  con¬ 
stant  at  19  Hz  with  the  amplitude  increasing.  At  this  time  the  aircr^t  was 
approaching  from  the  far  field:  for  all  practical  purposes  from  the  acous¬ 
tic  sensor’s  point  of  view,  the  aircraft  was  heading  straight  toward  the  sen¬ 
sor.  The  next  section  will  present  a  more  detailed  calculation  of  this 
frequency. 

From  65  sec  to  80  sec,  the  frequency  decreases  steadily  while  the 
amplitude  rises  to  its  maximum  and  then  starts  to  fall.  In  this  period  (at 
approximately  71  sec)  the  helicopter  passes  CPA.  The  slow  drop  in  fre¬ 
quency  guarantees  that  the  aircraft  will  pass  the  sensor  at  some  non-zero 
range.  The  slower  the  drop  in  frequency,  the  farther  the  range  will  be  at 
CPA.  The  key  piece  of  information  necessary  in  this  segment  is  the  slope 
of  the  frequency  curve. 

From  80  sec  to  120  sec,  the  frequency  remained  approximately  con¬ 
stant  again  at  between  15  and  16  Hz.  The  amplitude  was  decreasing  also. 
At  this  time  the  aircraft  was  in  the  far-Held  retreating  configuration:  from 
the  acoustic  sensor’s  point  of  view,  the  aircraft  was  going  directly  away 
from  the  sensor.  Since  the  frequency  falls  almost  halfway  between  two  in¬ 
teger  frequency  values,  the  accurate  calculation  of  this  frequency  can 
result  in  problems  concerning  which  frequency  bin  to  select. 


Approaching  Fundamental  Signal  - 19.10  Hz 

When  one  considers  the  approaching  signal  in  the  time  interval, 

[30,45],  from  30  to  45  sec  into  Run  21,  the  signal-to-noise  ratio  is  high, 
and  the  frequency  still  appears  to  be  constant.  This  is  the  appropriate  time 
to  determine  the  far-field  approaching  frequency. 


Complex  plane  of  the  FFT 

Table  5  shows  the  FFT  amplitudes  near  the  end  and  coherence  over  the 
entire  time  interval  for  a  number  of  frequencies  close  to  the  dominant  en¬ 
ergy  at  19  Hz.  Although  the  two  largest  amplitudes  are  usually  the  19  and 


The  units  of  the  amplitude  are  a  scalar  multiple  of  microbars.  The  scalar  was  applied 
to  obtain  the  maximum  dynamic  range  of  signal  from  the  analog  to  digital  processing  step. 
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Table  5 

Far-Field  Approaching  FFT  Amplitude  of  a  Helicopter 

Frequency 

(Hi) 

Coherence 

130,45] 

Time  Intervals  (see)  of  FFT  Analysis  | 

[38,30] 

[30,40] 

[40,41] 

[41,42] 

[42,43] 

[43,44] 

10.0 

0.178 

0.007 

0.012 

0.014  * 

0.022 

0.007 

0.029 

11.0 

0.002 

0.044 

0.019 

0.030 

0.011 

0.010 

0.015 

12.0 

0.002 

0.046 

0.014 

0.014 

0.012 

0.008 

0.011 

13.0 

0.047 

0.025 

0.021 

0.010 

0.009 

0.004 

0.027 

14.0 

0.040 

0.019 

0.005 

0.009 

0.002 

0.011 

0.027 

15.0 

0.240 

0.011 

0.016 

0.013 

0.004 

0.013 

0.002 

16.0 

0.291 

0.012 

0.018 

0.004 

0.011 

0.013 

0.010 

17.0 

0.088 

0.018 

0.017 

0.004 

0.012 

0.014 

0.014 

18.0 

0.047 

0.015 

0.009 

0.012 

0.008 

0.016 

0.010 

19.0 

0.700 

0.071 

0.082 

0.100 

0.105 

0.106 

0.070 

20.0 

0.160 

0.018 

0.008 

0.017 

0.017 

0.035 

0.008 

21.0 

0.080 

0.014 

0.004 

0.010 

0.011 

0.010 

0.012 

22.0 

0.175 

0.003 

0.006 

0.002 

0.003 

0.002 

0.003 

23.0 

0.000 

0.006 

0.001 

0.005 

0.002 

0.005 

0.001 

24.0 

0.032 

0.011 

0.006 

0.002 

0.005 

0.005 

0.007 

25.0 

0.011 

0.012 

0.005 

0.006 

0.004 

0.007 

0.001 

26.0 

0.083 

0.004 

0.001 

0.003 

0.004 

0.006 

0.002 

27.0 

0.056 

0.000 

0.003 

0.003 

0.001 

0.004 

0.003 

26.0 

0.126 

0.006 

0.005 

0.005 

0.009 

0.002 

0.005 

29.0 

0.171 

0.008 

0.003 

0.005 

0.004 

0.002 

0.005 

30.0 

0.321 

0.007 

0.007 

0.004 

0.001 

0.006 

0.009 

20  Hz,  the  coherence  .700  of  the  19-Hz  energy  is  much  higher  than  the 
20  Hz  coherence  of  .160.  (It  must  be  noted  that  the  coherence  was  calcu¬ 
lated  over  all  fifteen  1-sec  intervals  from  30  to  45  sec.  Hence,  it  appears 
to  contain  much  more  information  than  shown  in  Table  5.)  The  phase 
(Table  6)  appears  to  be  much  more  challenging  to  interpret  (even  for  the 
strong  19-Hz  energy).  The  plot  of  the  FFT  of  the  19-Hz  energy  in  the 
complex  plane  is  shown  in  Figure  12  for  the  time  period  from  30  to 
45  sec.  This  graph  shows  the  rotation  of  the  successive  points  around  the 
origin  in  a  manner  reminiscent  of  a  pure  cosine  wave  (see  Figure  7). 
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Table  6 

Far-Fleld  Approaching  FFT  Phase  (degrees)  of  a  Helicopter 

Frequency 

(Hz) 

Coherence 

130,45] 

Time  Intervals  (sec)  of  FFT  Analysis 

[38,39] 

[39,40] 

[40,41] 

[41,42] 

[42,43] 

[43,44] 

14.0 

0.040 

68.7 

129.5 

176.4 

-53.6 

-90.0 

13.1 

15.0 

0.240 

-113.4 

-131.2 

-168.2 

111.4 

-95.8 

128.8 

16.0 

0.291 

-143.9 

-95.2 

19.8 

-9.3 

-58.1 

110.7 

17.0 

0.088 

-145.8 

62.8 

-47.7 

134.0 

26.4 

3.6 

18.0 

0.047 

-79.2 

146.4 

-135.0 

-114.3 

-75.2 

-175.3 

19.0 

0.700 

-154.7 

-139.4 

-113.1 

-89.2 

-57.0 

-20.9 

20.0 

0.160 

-133.0 

-176.8 

117.5 

104.9 

113.7 

39.1 

21.0 

0.080 

-41.7 

89.4 

91.7 

92.4 

117.5 

-1.5 

22.0 

0.175 

-0.7 

-11.9 

5.0 

65.1 

90.8 

-8.5 

23.0 

0.000 

10.8 

-79.9 

53.5 

-41.5 

140.0 

-94.7 

24.0 

0.032 

37.9 

-69.7 

4.4 

155.5 

139.4 

125.6 

25.0 

0.011 

21.1 

3.3 

145.1 

-55.5 

147.7 

39.3 

Complex  plane  of  the  CSD 

The  CSD  of  successive  time  intervals  is  plotted  in  Figure  13.  The  CSD 
is  obtained  by  multiplying  the  two  FFT  amplitudes  and  taking  the  differ¬ 
ence  between  the  two  FFT  phases  of  adjacent  time  intervals.  The  roughly 
constant  rotation  of  the  FFT  phase  is  translated  to  the  CSD  as  the  points 
roughly  all  plotting  at  a  constant  angle  (the  rotation  angle)  to  the  positive 
real  axis.  All  of  the  large  amplitude  points  fall  in  the  first  quadrant  (posi¬ 
tive  real  and  imaginary  parts).  This  fact  is  exploited  when  a  number  of  ad¬ 
jacent  CSD  values  are  added  together. 


Cumulative  sum  of  the  CSD 

Figure  14  plots  the  cumulative  sums  of  the  CSD  values  of  adjacent  in¬ 
tervals:  The  first  point,  the  CSD  value  that  compares  [0,1]  to  [1,2],  is  con¬ 
nected  to  the  point  that  is  the  complex  sum  of  the  first  two  CSD  values. 
The  line  that  connects  these  two  points  actually  represents  the  second 
CSD  complex  number.  This  process  is  continued  until  all  14  of  the  CSD 
values  are  added  together.  Table  7  summarizes  the  calculation  of  the 
cumulative  sum  of  the  CSD  for  the  19-Hz  energy. 
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Table  7 

Calculation  of  the  Cumulative  Sum  of  the  Cross-Spectral 

Density  for  19  Hz  Far-Fieid  Approaching  Helicopter 

Time 

Intervals 

Cross-Spectral  Density 

Cumulativs  Sum 

First 

Second 

Real 

Imaginary 

Rsal 

Imaginary 

31 

[30.31] 

[31.32] 

0.000435 

-0.000545 

0.000435 

-0.000545 

32 

[31.32] 

[32.33] 

-0.000395 

0.001141 

0.000040 

0.000596 

33 

[32.33] 

-  1 

[33.34] 

0.001673 

0.003371 

0.001713 

0.003966 

34 

[33.34] 

[34.35] 

0.005093 

0.000698 

0.006806 

0.004664 

35 

[34.35] 

[35,36] 

0.003188 

0.005995 

0.009994 

0.01 06% 

36 

[35.36] 

[36,37] 

0.001883 

-0.001213 

0.011877 

0.009446 

37 

[36.37] 

[37.38] 

0.000653 

0.002739 

0.012530 

0.012185 

38 

[37.38] 

[38,39] 

0.003143 

0.005400 

0.015673 

0.017585 

39 

[38.39] 

[39,40] 

0.005616 

0.001536 

0.021288 

0.019121 

40 

[39.40] 

[40.41] 

0.007351 

0.003633 

0.026639 

0.022755 

41 

[40.41] 

[41.42] 

0.009600 

0.004254 

0.038239 

0.027009 

42 

[41.42] 

[42,43] 

0.009416 

0.005931 

0.047657 

0.032940 

43 

[42.43] 

[43.44] 

0.005995 

0.004372 

0.053652 

0.037311 

44 

[43.44] 

[44.45] 

0.002681 

0.004837 

0.056333 

0.042148 

Interpretation  of  approaching  sound  frequency 

The  final  point  X=.056333  +  .042148/  of  Figure  14  represents  the  vec¬ 
tor  sum  of  the  real  and  imaginary  components  of  the  14  individual  CSD 
values.  The  direction  of  this  point  from  the  origin  represents  the  CSD 
amplitude  weighted  sum  of  angles  of  the  individual  CSD  vectors  and  is, 
therefore,  the  best  representation  of  the  change  of  the  approaching  signal 
from  a  pure  19-Hz  component.  Using  Equation  20,  the  approaching  far- 
field  frequency  f,  can  be  calculated  more  accurately  as 


fa  =  19.00 //z  + 


=  19.00/fz  + 


f0.042l48'\ 

arctan 

0.056333 

_ _ __Z _ 

1  sec 

1  sec 


(23) 


=  19.10 //z 
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The  value  =  19.10  Hz  can  be  considered  as  the  average  frequency  of 
the  approaching  coherent  energy  falling  near  19  Hz  over  the  time  period 
[30,45].  The  degree  of  accuracy  is  presently  undetermined,  but  the  good 
alignment  of  the  individual  vectors  that  form  the  cumulative  sum  (Fig¬ 
ure  14)  suggests  that  the  above  number  is  accurate  to  within  0.1  Hz.  (The 
additional  precision  implied  in  19.10  Hz  will  be  retained  because  of  the 
sensitive  nature  of  some  future  calculations.) 

Another  test  of  accuracy  uses  Equation  20  on  each  of  the  14  CSD  meas¬ 
urements  to  determine  how  much  the  actual  frequency  implied  by  the 
phase  of  the  individual  CSD  values  varies  over  that  length  of  time.  Fig¬ 
ure  IS  summarizes  the  results.  Most  of  the  individual  frequencies  fall  be¬ 
tween  19.0  and  19.2  Hz.  The  Hrst  two  CSD  values  (when  20  Hz  had  a 
larger  amplitude)  increased  the  mean  to  19.15  Hz  and  the  standard  devia¬ 
tion  to  .16  Hz.  (Ignoring  those  first  two  CSD  values,  the  mean  and  stand¬ 
ard  deviation  become  19.10  and  .08  Hz,  respectively.  This  would  agree 
with  Equation  23,  which  is  a  19-Hz-amplitude  weighted  sum  of  CSD  com¬ 
plex  numbers.) 


Comparison  with  neighboring  frequencies 

Amplitudes.  The  main  energy  of  the  CSD  and  FFT  is  concentrated  at 
19  Hz  in  the  time  period  [33,45]  (Table  5).  The  fact  that  the  19-Hz  FFT 
amplitude  is  over  three  times  larger  than  the  amplitudes  of  15  to  18  and  20 
to  30  Hz  in  [38,44]  suggests  that  analyzing  the  19-Hz  and  CSD  infor¬ 
mation  is  sufficient  to  determine  the  approaching  frequency  of  the  source 
(Figure  16). 

Coherences.  The  coherence  (Table  5  and  Figure  17)  of  the  1 -sec -long 
intervals  composing  the  total  interval  [30,45]  is  another  good  indicator 
that  the  19-Hz  frequency  is  the  only  FFT  and  CSD  energy  that  should  be 
considered.  Of  all  the  frequencies  between  10  and  30  Hz,  only  the  19-Hz 
component  had  a  coherence  above  0.325  over  the  entire  15-sec-long  time 
peri^  [30,45].  At  0.700,  the  19-Hz  coherence  was  more  than  twice  as 
large  as  the  maximum  of  those  other  frequencies.  (Because  the  18-Hz 
coherence  is  0.047  and  the  20  Hz  is  0.160,  far  below  the  19-Hz  level, 
these  frequencies  must  not  be  used  in  the  far-field  analysis.  This  result 
for  the  approaching  signal  was  significantly  different  from  the  situation 
for  the  retreating  signal,  where  frequencies  adjacent  to  the  maximum 
amplitude  frequency  also  needed  to  be  analyzed.) 


Retreating  Fundamental  Signal  •  15.65  Hz 

The  retreating  far-field  signal  vastly  differed  from  the  approaching  far- 
field  signal.  As  Figure  11  indicates,  the  approaching  signal  was  evident  at 
about  60  sec  before  the  CPA.  However,  the  retreating  signal  has  dropped 
off  the  threshold  of  the  graph  within  30  sec  after  the  CPA.  For  this  reason. 
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the  analysis  of  the  retreating  signal  was  performed  much  closer  to  the  CPA 
in  the  interval  [90,105]. 


Comparison  of  neighboring  frequencies 

Ampiitudes.  The  strong  energy  in  this  time  interval  was  around  IS 
and  16  Hz.  Tlie  maximum  energy  was  at  16  Hz,  with  15  Hz  in  second 
place  at  about  50  percent  of  maximum;  third  was  17  Hz  at  roughly  30  per¬ 
cent  of  maximum.  Table  8  summarizes  the  amplitude  for  a  selected  por¬ 
tion  of  this  time  period.  The  corresponding  phase  information  is  given  in 
Table  9. 


Tabie  8 
Far-Field 

Retreating  FFT  Amplitude  of  a  Moving  Helicopter 

Frequency 

(Hz) 

Coherence 

[90,1051 

Time  Intervale  (sac)  of  FFT  Analyala 

[90,911 

[91,921 

[92,931 

[93,941 

[95,961 

[97,981 

10.0 

0.150 

0.013 

0.015 

0.013 

0.008 

0.004 

0.006 

11.0 

0.051 

0.004 

0.007 

0.019 

0.009 

0.005 

0.013 

12.0 

0.088 

0.010 

0.010 

0.008 

0.011 

0.003 

0.012 

13.0 

0.632 

0.011 

0.013 

0.016 

0.010 

0.011 

0.009 

14.0 

0.607 

0.017 

0.014 

0.027 

0.019 

0.020 

0.021 

15.0 

0.717 

0.056 

0.054 

0.040 

0.051 

0.041 

0.055 

16.0 

0.877 

0.109 

0.101 

0.109 

0.067 

0.081 

0.082 

17.0 

0,712 

0.034 

0.034 

0.023 

0.020 

0.025 

0.028 

16.0 

0,376 

0.021 

0.008 

0.011 

0.012 

0.014 

0.027 

19.0 

0.371 

0.047 

0.027 

0.026 

0.007 

0.032 

0.027 

20.0 

0.170 

0.041 

0.029 

0.018 

0.028 

0.033 

0.010 

21.0 

0.429 

0.013 

0.015 

0.011 

0.005 

0.010 

0.003 

22.0 

0.402 

0.009 

0.013 

0.004 

0.011 

0.005 

0.004 

23.0 

0.150 

0.006 

0.007 

0.005 

0.011 

0.005 

0.007 

24.0 

0.179 

0.004 

0.010 

0.006 

0.005 

0.001 

0.006 

25.0 

0.002 

0.002 

0.012 

0.006 

0.006 

0.001 

0.007 
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Table  9 

Retreating  FFT  Phase  (degrees)'  of  a  Moving  Helicopter 

Frequency 

(Hz) 

Coherence 

(00,105] 

Time  Intervals  (sec)  of  FFT  Analysis 

[90,91] 

[91.92] 

[02.03] 

[93.94] 

[95,96] 

[97,98] 

10.0 

0.179 

•116.4 

102.8 

30.2 

-36.0 

157.7 

-21.7 

11.0 

0.094 

-98.9 

152.3 

-1.4 

-15.2 

-173.1 

86.5 

12.0 

0.141 

-123.2 

164.9 

4.1 

-17.3 

131.9 

45.4 

13.0 

0.710 

-120.6 

119.7 

40.6 

-69.6 

147.5 

51.9 

14.0 

0.709 

-101 .6 

153.7 

10.7 

-71.1 

149.3 

48.1 

15.0 

0.768 

-114.9 

144.5 

3.2 

•82.1 

176.9 

40.9 

16.0 

0.893 

71.1 

-50.0 

•152.9 

92.2 

-11.6 

-144.2 

17.0 

0.781 

72.0 

-52.9 

-137.1 

132.6 

-13.7 

-109.3 

18.0 

0.593 

46.2 

-40.1 

-154.9 

72.0 

0.1 

-129.4 

19.0 

0.580 

47.0 

-161.3 

-85.1 

•26.0 

80.6 

-112.7 

20.0 

0.261 

165.4 

20.7 

155.9 

53.5 

•64.3 

56.1 

21.0 

0.581 

130.7 

-6.4 

-178.6 

76.8 

-27.8 

165.1 

22.0 

0.644 

153.3 

•23.2 

•175.1 

68.9 

-39.6 

130.2 

23.0 

0.186 

-143.0 

12.6 

-179.8 

87.9 

•64.1 

171.7 

24.0 

0.254 

82.4 

-33.7 

-178.9 

122.7 

88.3 

-177.9 

25.0 

0.003 

-79.6 

-48.1 

•169.3 

-24.6 

-11.5 

-170.2 

'  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI  units  is  presented  on 
page  vi. 

Coherences.  Unlike  the  approaching  signal,  the  coherence  remained 
above  0.70  for  all  three  of  the  FFT  frequencies  of  15,  16,  and  17. 

Cumulative  sums  of  the  CSD.  Because  of  the  adequate  amplitude  and 
good  coherence,  all  three  of  these  frequencies  were  analyzed  using  the 
cumulative  sums  of  CSD  (Figure  18).  If  the  different  amplitudes  are  ig¬ 
nored,  all  three  frequencies  show  approximately  the  same  shaped  curve. 


Frequency  calculations 

Table  10  lists  the  actual  cumulative  sum  of  CSD  of  individual  FFT  fre¬ 
quencies,  coherence,  and  the  implied  frequency  change  calculated  from 
l^uation  20. 


Chaptsr  3  Analysis  of  Far-Rekl  Frequency  for  Measured  Data 


27 


Table  10 

Calculation  of  Retreating  Far-Field  Frequency 

Frequency 

hertx 

Coherence 

Cumulative  Sum  CSD 

Implied  Angle  Change 

Real 

Imaginary 

degrees 

cycles 

15 

0.717 

-0.011426 

-0.013759 

-129.71 

-0.360 

16 

0.877 

-0.032258 

-0.049309 

-123.19 

-0.342 

17 

0.712 

-0.003260 

-0.005241 

-121.88 

-0.339 

Vector  Sum 

-0.046944 

-0.068309 

-124.50 

-0.346 

interpretation 

Since  all  of  the  three  significant  frequencies  look  similar,  the  vector 
sum  (Figure  19)  probably  represents  the  best  approximation  of  the  total 
signal.  Note  that  this  sum  applies  weights  appropriate  to  the  respective 
amplitudes  of  each  of  the  individual  signals.  Since  16  Hz  has  the  largest 
amplitude,  most  likely  the  fmal  change  number  represents  the  shift  from 
16  Hz. 

/r  =  retreating  frequency  =  16  -  0.35  =  15.65^z  (24) 

Because  the  absolute  value  of  the  angle  change  is  close  to  0.5  cycle,  it 
is  best  to  check  the  frequency  of  the  second-largest  amplitude  (to  guaran¬ 
tee  that  the  answer  is  not  16.65  Hz).  Since  the  15-Hz  amplitudes  are  sec¬ 
ond  largest  behind  the  16  Hz,  the  energy  can  be  expected  to  fall  between 
15  and  16  Hz  as  the  calculated  fj  satisfies. 
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4  Combination  of  Both 
Frequencies  for  the 
Fundamental 


Basic  Frequency  Calculation 

The  combined  knowledge  of  both  the  approaching  and  retreating  fre¬ 
quencies  permits  the  calculation  of  both  the  basic  frequency  of  the  helicop¬ 
ter  and  its  velocity  (under  the  assumptions  of  Table  1).  The  values  of  f^ 
and  ff  (Equations  23  >>'  J '  >)  can  be  inserted  into  Equation  3  to  yield  an 
estimate  of  the  basic  fr;  .^uency  fg  that  the  helicopter  would  generate  if  it 
were  not  moving: 


/  U±Il  2  •  19.10-  15.65 
*  19.10+  15.65 


1 

=  17.21  Hz 


(25) 


This  corresponds  to  the  reported  value  of  17.2  Hz,  four  times  the  fre¬ 
quency  of  rotation  of  the  main  rotor  Shaft.  With  four  blades  on  that  rotor, 
the  dominant  energy  in  Figure  1 1  appears  to  be  caused  by  the  main  rotor 
blades.  The  passive  acoustic  sensor  analysis  has  correctly  detected  this 
frequency. 


Since  the  microphone  is  measuring  this  frequency  when  the  helicopter 
is  at  CPA,  Figure  5  indicates  that  CPA  occurred  sometime  around  70  sec 
into  Run  2 1 . 
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Velocity  Estimation 

Likewise,  the  velocity  of  the  source  can  be  estimated  using  Equation  4 
or  6. 


V  =  C  ■ 


fa-fo 

fa 


=  352.6  • 


19.10-  17.21 
19.10 


(26) 


=  35.0 

The  formula  requires  only  the  knowledge  of  the  approaching  frequency 
fg,  the  basic  frequency  f^,  and  the  velocity  of  sound  C.  This  formula  im¬ 
plies  that  if  the  basic  frequency  f^  were  known  while  the  helicopter  was 
still  in  the  far-field-approaching  situation,  the  velocity  could  be  calculated 
as  soon  as  the  approaching  far-Held  frequency  f,  could  be  determined.  In 
Run  21  this  frequency  was  calculated  using  the  data  from  30  to  45  sec 
(see  "Approaching  Fundamental  Signal  -  19,10  Hz,"  pp  21-25)  long  before 
the  helicopter  reached  CPA  at  after  70  sec.  Hence  the  velocity  of  the  in¬ 
coming  helicopter  could  be  estimated  over  25  sec  before  it  reached  CPA. 
The  actual  number  compares  favorably  with  the  velocity  of  65  knots 
(33.5  m/sec)  reported  by  the  pilots.  These  results  apply  to  this  data  set. 
The  accuracy  of  the  approaching  velocity  estimate,  the  time  before  CPA 
that  it  can  be  calculated,  and  the  length  of  the  time  intervals  required  in 
the  calculation  are  all  dependent  on  the  amount  of  background  noise,  the 
amplitude  of  the  incoming  signal,  the  quality  of  acoustic  recording  equip¬ 
ment,  and  the  possible  presence  of  other  targets.  The  u.  /estigation  of  the 
limitations  of  the  velocity  formula  is  a  fertile  area  for  study  — but  remains 
outside  the  scope  of  this  report. 
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5  Analysis  of  the  CPA 
Frequencies  for  the 
Fundamental 


Calculation  of  the  Frequency  at  One  Time 

Witii  the  estimates  of  basic  frequency  f^  ^  17.21  Hz  from  Equation  25 
and  the  velocity  of  the  helicopter  V  =  35.0  m/sec  (Equation  26),  the  time 
of  CPA  and  range  at  CPA  can  be  estimated  from  the  CSD  data  of  Table  11. 
The  aircraft  was  at  CPA  when  it  generated  a  frequency  f^  at  the  receiver. 

The  FFT  Amplitude  part  of  Table  1 1  suggests  that  the  signal  crosses 
through  17.21  Hz  during  the  interval  [71,72J.  TTierefore,  the  detailed 
analysis  of  the  frequency  using  Equation  20  on  the  17-Hz  component  of 
the  CSD  that  compares  the  [70,71]  to  [71,72]  should  generate  the  range  es¬ 
timate.  The  results  will  be  considered  applicable  to  71  sec.  For  17  Hz: 

CSD  Ampl  71  =  (FFT  Ampl  [70,71])  (FFT  Ampl  [71,72]) 

=  (0.407)  (0.873)  =  0.355311 

CSD  Phase  71  =  (FFT  Phase  [71,72])  -  (FFT  Phase  [70,71]) 

=  (-104.9)  -  171.3  deg  =  -276.2  deg 
=  360.0  -  276.2  deg  =  83.8  deg 
=  83.8  deg  (1  cycle/360  deg) 

=  0.23  cycle 

Since  the  17-Hz  CSD  amplitude  is  an  obvious  maximum  for  all  of  the 
CSD  frequencies  from  15  Hz  to  20  Hz  of  the  time  71  sec.  Equation  20  sug¬ 
gests  that  the  dominant  frequency  f(71)  of  the  signal  at  71  sec  is 
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y(71)  =  17  //2  +  =  17.23  Hz 

1  sec 


The  value,  calculated  from  the  CSD  values  of  [70,71]:[71,72],  can  be 
checked  against  the  original  RFT  results  for  [70,71]  and  [71,72].  The 
amplitude  of  the  FFT  of  [70,71]  is  .407  at  17  Hz  and  .339  at  18  Hz  (Table 
11),  suggesting  that  the  main  frequency  at  time  70.5  is  roughly  halfway  be¬ 
tween  17  and  18  Hz  —  but  closer  to  17  Hz.  So  the  FFT  of  [70,71]  sug¬ 
gests  that  f(70.5)  e  17.4  Hz.  Similarly,  the  FFT  of  [71,72]  suggests  that 
f(7l.S)  B  17  1  Hz,  since  the  amplitude  of  18  Hz  is  slightly  larger  than  that 
of  16  Hz.  Taking  the  average  of  these  two  approximations  of  the  frequen¬ 
cies  of  the  FFT’s  yields  an  average  of  17.25  Hz  for  an  f(70).  Hence  the  in¬ 
dividual  FFT  amplitudes  of  the  1 -sec-long  intervals  compare  well  with  the 
CSD  calculation  of  17.23  Hz  for  the  boundary  between  the  time  intervals. 


Frequency  Estimate  as  a  Function  of  Time 


The  interpretation  of  the  frequency  of  the  dominant  energy  was  made 
for  all  of  the  times  from  61  to  74  sec  (see  Figure  20).  The  calculation  of 
the  frequencies  from  68  to  72  sec  is  shown  in  Table  11.  These  calcula¬ 
tions  were  made  using  the  FFT  frequency  at  the  largest  amplitude.  The 
smooth  graph  in  Figure  20  suggests  that  the  frequency  calculation  is 
robust  (i.e.,  not  subject  to  large  errors  resulting  from  noise  in  the  data). 


The  time  t^  when  the  sound  reaches  the  receiver  of  the  helicopter  at 
CPA  can  be  estimated  using  Equation  10  and  the  straight  lines  through  suc¬ 
cessive  points  t-1  and  t  for  the  time  increment  of  1  sec  of  Figure  20.  This 
calculation  is  based  on  the  intersection  of  the  basic  frequency  with  the 
linear  interpolation  of  the  Doppler-shifted  frequency  as  a  function  of  time. 


Ut)  =  t  -I- 


fo-m 


dt 


(0 


toit)  =  t  + 


17.21  Hz  -  At) 

m  -  -  1) 


•  (1  sec) 


This  formula  is  evaluated  in  Table  12  using  the  interpreted  frequencies 
from  the  bottom  of  Table  11.  All  of  the  values  are  close.  The  time  of  CPA 
of  71.07  sec  will  be  used  for  the  estimate  for  the  main  rotor  blades.  Since 
this  time  was  reasonably  estimated  before  the  helicopter  reached  CPA,  the 
time  of  CPA  can  be  calculated  sufficiently  early  to  decide  on  a  response  to 
this  incoming  target. 
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Table  12 

Estimates  of  the  Time  of  CPA 

Time  t  of  Second  Point 
of  Slope,  sec 

Slope,  Hx/sec 

Estimate  of  CPA  to(t),  sec 

69 

-0.38 

70.92 

70 

-0.34 

71.11 

71 

-0.35 

71.07 

72 

-0.36 

71.07 

Calculation  of  the  Range  Using  the  Slope 

The  slope  of  the  frequency  as  a  function  of  time  (Figure  20)  appears  to 
remain  roughly  constant  from  68  sec  to  74  sec.  This  means  that  the  for¬ 
mula  for  the  range  at  CPA,  which  requires  the  slope  at  CPA,  can  be 
evaluated  using  the  approximations  of  the  slope  before  the  helicopter 
reaches  CPA.  The  main  rotor  prediction  line  of  Figure  2 1  is  the  applica¬ 
tion  of  Equation  5 


C  ■  ^  iCPA  ) 

with  current  estimates  of  the  data 


R(t)  =  -17.21  Hz  •  (35.0  ysec?  (27) 

352.6  '%ec  •  (/(f)  -  /(f  -  1)1  /  (1  sec) 

The  values  near  CPA,  summarized  in  Table  13,  are  stable  in  the  vicinity 
of  167  m  even  before  the  helicopter  arrives  at  CPA.  The  estimated  range 
given  by  the  pilots  was  182.3  m  [100  m  east  (estimated)  and  152.4  m 
vertical]. 


Summary  of  the  Analysis  of  the  Main  Rotor 
Signal 


The  CSD  analysis  of  the  strong  main  rotor  energy  that  initially  fell 
around  19  Hz  yielded  a  great  deal  of  relevant  information  about  the  heli¬ 
copter.  The  cumulative  sum  of  the  cross-spectral  densities  over  the  time 
intervals  [30,45]  and  [90,105]  yielded  estimates  of  the  approaching  and 
retreating  frequencies  at  19.10  and  15.65  Hz,  respectively.  These  values 
yielded  the  estimate  of  the  basic  frequency  of  17.21  Hz,  which  agrees 
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Table  13 

Estimates  of  the  Range  at  CPA 

Time  t  of  Second 
Point  of  Slope,  sec 

Estimate  of  CPA 
Time,  sec 

Estimate  of  Range,  m 

Error  of  Range 

69 

70.92 

158.4 

-12.9% 

70 

71.11 

174.6 

-4.1% 

71 

71.07 

168.7 

-7.5% 

72 

71.07 

165.0 

-9.4% 

with  the  known  17.2-Hz  frequency  of  the  four  rotor  blades  of  the  helicop¬ 
ter.  The  basic  frequency  and  the  approaching  frequency  were  combined  to 
calculate  the  velocity  of  the  aircraft  as  35.0  tn/sec,  slightly  above  the 
reported  33.5  m/sec  reported  by  the  pilots  (but  less  than  some  preliminary 
global  positioning  system  (GPS)  data). 

Using  the  basic  frequency  and  the  slope  of  the  Doppler-shifted  fre¬ 
quency  as  a  function  of  time  curve,  the  time  of  CPA  is  estimated  to  be  at 
approximately  71  sec  into  the  run  (Figure  20).  The  basic  frequency,  slope 
near  CPA,  and  velocity  estimates  combine  to  predict  the  range  167  m  at 
CPA  (Figure  21).  The  rough  estimate  of  the  pilots  was  182  m. 

The  CSD  analysis  estimated  within  acceptable  tolerances  the  basic  fre¬ 
quency  (of  the  fundamental),  velocity,  time  of  CPA,  and  range  of  CPA  of 
the  moving  helicopter.  It  is  very  significant  that  all  of  this  information 
was  derived  solely  from  passive  acoustic  measurements  at  one 
microphone. 
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6  Conclusions  and 
Recommendations 


Conclusions 

Passive  acoustic  energy  can  be  used  to  estimate  the  velocity,  range  at 
CPA,  and  time  of  CPA  from  knowledge  of  the  Doppler  shift  of  frequency 
over  time  recorded  at  one  acoustic  receiver  (Equations  S,  6,  and  11).  For¬ 
mulas  for  these  estimates  require  a  high  resolution  of  the  frequency  of  the 
signal  over  relatively  short  time  periods. 

The  traditional  FFT  frequency  analysis  is  not  sufficient  in  this  dynamic 
situation;  the  CSD  is  shown  to  be  much  better  than  conventional  FIT  ap¬ 
proaches.  The  phase,  measured  in  cycles  and  divided  by  the  length  of  the 
FFT  interval,  can  be  interpreted  as  the  amount  of  the  shift  between  the  ac¬ 
tual  dominant  energy  and  the  frequency  of  the  FFT  calculation  (Equa¬ 
tion  20).  For  the  dominant  signal  of  the  data,  the  CSD  frequency 
resolution  appeared  to  be  within  0.1  Hz. 

When  the  basic  frequency  (i.e.,  frequency  before  any  Doppler  shift)  is 
known,  the  high  resolution  of  the  far-Held  approaching  frequency  of  the 
dominant  main  rotor  energy  at  19  Hz  permits  a  prediction  of  35.0  m/sec 
for  the  approach  velocity,  accurate  to  within  S  percent  of  the  value  of 
33.5  m/sec  reported  by  the  pilots.  The  analysis  of  this  energy  near  CPA 
yielded  range  estimates  of  158  to  175  m,  accurate  to  within  13  percent  of 
the  range  of  182  m  reported  by  the  pilots.  Some  of  the  sensor  estimates 
were  made  before  the  helicopter  reached  CPA. 

If  nothing  were  known  about  the  approaching  aircraft,  the  basic  fre¬ 
quency  and  all  of  the  above  information  could  be  calculated  to  the  same 
levels  of  accuracy  by  analyzing  the  far-Held  retreating  frequency  after  the 
aircraft  had  moved  far  away.  The  combined  knowledge  of  the  far-field  ap¬ 
proaching  frequency  and  the  far-field  retreating  frequency  yields  an  es¬ 
timated  basic  frequency  of  17.21  Hz.  The  specification  of  the  helicopter 
is  17.2  Hz.  The  presence  of  two  sensors,  with  one  far  out  in  front  of  the 
other,  can  provide  resolution  of  the  unknown  frequency  problem.  If  the 
sensor  out  front  has  the  task  of  determining  the  basic  frequency,  then  this 
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frequency  can  be  passed  to  the  second  sensor  to  determine  the  range  when 
the  aircr^t  approached  CPA. 

Calculations  provided  herein  used  the  data  from  one  passive  sensor  on 
actual  acoustic  data.  The  first  harmonic  and  the  tail  rotor  energy  provided 
additional  measurements  that  represent  independent  statistical  support  of 
the  main  rotor  calculations  of  velocity  and  range.  Additional  sensors  add 
independent  data.  All  of  these  measurements  can  greatly  enhance  the 
capability  of  using  acoustic  data  to  perform  as  a  passive  detection  and 
ranging  device. 


Recommendations 

The  following  is  a  list  of  topics  recommended  for  future  study: 

a.  Creation  of  a  computer  system  to 

(1)  detect  and  track  coherent  energy. 

(2)  perform  the  above  analysis  automatically. 

(3)  perform  calculations  on  an  80386  PC  in  real  time. 

b.  Integration  of  the  computer  system  with 

( 1 )  direction-fmding  programs . 

(2)  wide-area  mine  (WAM)  technology. 

c.  Further  study  of 

(1)  FFT  intervals  of  0.5  sec  and  0.25  sec. 

(2)  faster  moving  sources. 

{a)  close  to  the  speed  of  sound. 

(b)  faster  than  the  speed  of  sound. 

(3)  lower  flying  targets. 

(4)  changing  velocities  of  the  aircraft. 

(5)  impact  of  curvilinear  flight  paths. 

(6)  effect  of  wind. 

(7)  effect  of  multiple  sources. 
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(8)  impact  of  battlefield  noise  interference. 

d.  Identification  of  acoustic  sources  that  emit  a  temporally  coherent 
signal. 

e.  Creation  of  a  list  of  the  basic  frequencies  of  each  aircraft. 

Such  studies  could  lead  to  broad  applications  in  both  civil  and  military 
aircraft  identification  and  tracking,  and  enhance  current  procedures. 
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Figure  2.  The  key  input  data  for  passive  estimation  of  velocity,  range,  basic  frequency,  and  time  of  CPA  if  nothing  is  known  about  the 
helicopter 


Velocity  B  33.5  m/sec 


e  3.  Viewpoint  comparisons:  source  versus  receiver  from  simulation  of  Doppler  shift  -  Run  21 
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Figure  5.  Frequency  of  maximum  FFT  amplitude  (main  rotor  blades)  of  helicopter  -  Run  21 


Figure  7.  The  19-Hz  FFT  analysis  of  temporal  coherence  (input  Is  a  19.1-Hz  cosine  wave).  Note:  coherence  is  1.000.  The  label  for 
each  point  represents  the  end  of  the  1 -sec-long  interval  (e.g.,  7  represents  the  FFT  of  the  time  interval  [6,7]) 
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Figure  9.  The  10-Hz  cumulative  sum  of  CSD  (input  is  a  19.1 -Hz  cosine  wave).  Point  7  represents  the  sum  of  the  CSDs  from  [0.1] 
(1,2]  up  to  [6.7] :  [7,8].  The  coherence  is  0.881  for  this  frequency 


Vicksburg,  MS  Airport,  Moving  Helicopter  Test 


Road  Railroad  Tracks 


FFT  Amplitude  of  Moving  Helicopter 


Figure  1 1 .  FFT  analysis  of  a  moving  helicopter 
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Figure  12.  FFT  analysis  of  temporal  coherence  of  19-Hz  energy  of  approaching  helicopter  during  the  time  interval  1^,45].  The 
coherence  was  0.700  for  this  total  time  period  of  Run  21.  The  label  represents  the  ending  time  of  the  1 -sec-long  FFT  interval 


Figure  13.  Cross-spectral  density  analysis  (19  Hz,  approaching  helicopter.  Run  21).  The  point  labeled  42  is  the  CSD  comparison  of 
the  FFTs  of  the  interval  [41,42] :  [42,43]  and  is  an  estimate  of  the  19-Hz  signal  at  42  sec 
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Figure  15.  Frequency  predicted  by  individual  cross-spectral  densities  of  approaching  main  rotor  energy,  Run  21,  channel  1. 
Note:  the  standard  deviation  (StDev)  is  0.16  Hz 
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Figure  16.  Comparison  of  amplitudes  of  FFT  of  approaching  19-H2  energy.  Run  21,  channel  1  where  velocity  «  33  m/sec,  altitude 
-  152  m,  horizontai  offset  -  100  m  (perpendicular  offset  -  182  m).  The  time  is  the  end  of  the  1-sec-long  Interval 
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Figure  17.  Coherence  over  [30,45]  of  lar-field  approaching  energy,  Run  21,  channel  1  moving  helicopter 
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Figure  18.  Cumulative  sum  of  CSD  for  retreating  helicopter  of  frequencies  15, 16,  and  17  Hz  with  coherences  0.717,  0.877,  and  0.712 
over  time  [90,105]  -  Run  21,  Vicksburg  Airport,  June  21,  1990 
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Figure  21.  Prediction  of  range  at  CPA  utilizing  the  slope  of  Doppler  shift.  Run  21,  channel  1.  The  pilot-reported  range  includes  the 
0.52-sec  shift  from  CPA  to  receiver 


Appendix  A 
Definitions 


Term 

Description 

Alias 

The  process  by  which  a  continuous  anaiog  high  frequency  f 
(above  the  Nyquist  frequertcy  fn)  signal  appears  to  look  like  a 
lower  frequency  f,  after  digitally  sampling  at  discrete  times,  as 
determined  by  the  formula: 

t,  =  Minimum  \  t  ~  k  ■  Z  ■  f„\ 

where  the  minimum  is  computed  over  all  positive  integers  k. 

Basic  frequency 

Any  isolated  band  of  frequency  energy  that  has  not  been  Doppler 
shifted  (e.g.,  fundamental  of  main  rotor  blades,  second  harmonic 
of  main  rotor  blades,  third  harmonic  of  main  rotor  blades, 
fundamental  of  tail  rotor  blades, ...) 

Frequency  bin 

The  discrete  number  that  the  Fast  Fourier  Transform  (FFT) 
assigns  to  each  frequency  of  a  real  input  signal.  For  example,  a 

1 -sec-long  19.1 -Hz  signal  would  have  most  of  its  energy 
assigned  to  the  19-Hz  FFT  bin  with  the  second  largest  amount  of 
energy  assigned  to  the  20-Hz  FFT  bin. 

Fundamental  frequency 

The  lowest  frequency  of  a  periodic  function. 

Harmonic 

A  frequency  which  is  an  integer  multiple  of  a  fundamental 
frequency.  The  second  harmonic,  for  example,  has  a  frequency 
twice  that  of  the  fundamental. 

Nyquist  frequency 

The  maximum  frequency  fn  that  can  be  accurately  recorded 
digitally  at  a  particular  sampling  rate  dt  in  time: 

f  ~  ^ 

"  2  dt 

Periodic  function 

A  function  that  re(.  3ats  after  a  specific  interval  of  time. 

A  Definitions 


Appendix  B 
Clarification  of  Various 
Frequencies 


This  report  mentions  many  different  types  of  frequency  measurements 
and  variables.  The  key  to  understanding  is  to  realize  the  precise  type  of 
frequency  that  is  under  discussion  at  each  stage.  This  appendix  attempts 
to  clarify  these  concepts  as  they  relate  to  physics  (i.e.,  acoustics)  and  ap¬ 
plied  mathematics  and  computing  (i.e.,  signal  processing). 

Table  B1  separates  the  frequencies  into  four  categories.  The  first 
category  represents  the  mechanical  revolution  rate.  The  second  category 
applies  to  those  frequencies  that  can  be  acoustically  detected  at  the 
source.  This  category  includes  the  prominent  mechanical  frequencies  and 
harmonics.  Any  frequency  in  this  second  category  is  considered  a  “basic 
frequency”  in  this  document.  There  are  some  major  changes  in  the  fre¬ 
quencies  from  where  the  sound  is  generated  to  where  it  is  recorded  at 
some  stationary  receiver  on  the  ground.  The  third  category  describes  how 
the  acoustic  signal  is  sensed  at  the  microphone  on  the  ground.  Such  fre¬ 
quencies  are  usually  referred  as  “Doppler-  shifted  frequencies”  in  the  text. 
Tlie  Hnal  category  represents  a  leap  into  the  mathematical  world  of  the 
standard  computational  tool,  the  Fast  Fourier  Transform  (FFT),  that  is 
used  to  analyze  measured  signals  of  amplitude  as  a  function  of  time.  This 
world,  “the  frequency  domain,”  reports  ail  of  its  measurements  in  terms  of 
frequencies:  hence  it  takes  more  work  in  that  domain  to  detect  and 
analyze  the  acoustical  frequencies  that  can  be  measured  at  the  receiver. 


Mechanical  Frequencies 

The  mechanical  frequencies  are  those  that  are  most  easily  understood 
by  the  layman.  They  are  the  frequency  of  revolution  of  physical  moving 
parts  of  a  machine.  Hence,  they  can  be  seen  and  readily  comprehended. 
For  a  helicopter  (see  Table  B2)  they  include  the  revolution  rate  of  the 
rotor  shafts  and  engine.  The  blades  also  represent  a  revolution  rate  that  is 
the  product  of  the  shaft  revolution  rate  times  the  number  of  blades  on  that 
shaft.  This  formula  assumes  that  the  blades  are  all  designed  alike  and  are 
spaced  an  equal  distance  around  the  shaft. 
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Acoustically  Generated  Frequencies  at  Source 

Although  there  are  many  mechanical  frequencies  on  the  helicopter,  not 
all  of  them  generate  significant  acoustic  signal.  The  blades  generate  much 
more  noise  than  the  shafts.  For  the  purpose  of  this  report  the  shafts  do  not 
appear  to  generate  any  measurable  signal  and  will  not  be  considered  as  an 
acoustically  generated  frequency.  The  blade  frequencies  are  acoustically 
measurable,  however  (see  Table  B3). 

If  the  acoustical  signal  generated  by  one  set  of  blades  appears  to  be  a 
perfect  sine  wave,  then  reporting  that  frequency  is  sufficient.  In  most 
cases  the  wave  does  not  appear  to  be  a  perfect  sine  wave.  Higher  har¬ 
monics  (see  Appendix  A)  are  present  in  the  acoustical  signal.  Normally, 
each  higher  harmonic  has  a  lower  amplitude  than  its  predecessor. 

Any  acoustically  generated  frequency  at  the  source  is  synonymous  with 
the  term  “basic  frequency.”  As  documented  in  the  Doppler  theory  section 
of  (Chapter  2  of  this  report,  if  the  basic  frequency  is  known  and  not  chang¬ 
ing  over  time,  then  accurate  knowledge  of  the  Doppler-shifted  frequency 
recorded  at  a  microphone  Hxed  on  the  ground  as  Ae  helicopter  is  ap¬ 
proaching  permits  the  determination  of  the  velocity  and  range  at  the 
closest  point  of  approach  (CPA)  of  a  helicopter  traveling  at  a  constant 
velocity  in  a  straight  line.  If  the  basic  frequency  is  not  known,  it  can  be 
determined  through  the  knowledge  of  the  Doppler-shifted  frequency  at  the 
microphone  using  all  time  (both  approaching  and  retreating). 

The  basic  frequencies  are  measured  at  the  source  (also  referred  to  as 
“the  helicopter  point  of  view”).  They  should  also  be  sensed  by  a  micro¬ 
phone  that  is  close  to  the  helicopter  when  neither  helicopter  nor 
microphone  is  moving  and  there  is  no  wind. 


Acoustically  Measured  Frequencies  at  Receiver 

At  the  normal  receiver  (e.g.,  a  microphone)  the  recorded  signal  has 
been  changed  from  those  generated  at  the  source.  When  the  source  is 
moving  relative  to  the  receiver  and  there  is  no  wind,  a  Doppler  shift  in  the 
frequency  occurs: 


m  ^fo  • 


where 


(1) 


f(t)  =  Doppler-shifted  frequency  at  time  t 
fo  =  basic  frequency 

V(t)  =  velocity  that  the  source  is  approaching  the  receiver  at  time  t 


B2 


Appendix  B  Clarification  of  Various  Frequencies 


C  =  velocity  of  sound 

Figure  1  in  the  main  text  illustrates  how  the  Doppler-shifted  frequency 
compares  to  the  basic  frequency.  As  the  source  is  approaching  the 
receiver,  V(t)>0  and  the  Doppler-shifted  frequency  is  higher  than  the 
basic  frequency.  At  CPA,  V(t)=0  and  they  are  both  equal.  As  the  helicop¬ 
ter  is  going  away  from  the  microphone,  the  Doppler-shifted  frequency  is 
less  than  the  basic  frequency.  The  microphone  senses  the  Doppler-shifted 
frequency,  not  the  basic  frequency. 

Some  additional  changes  are  seen  in  the  recorded  signal:  first,  the 
amplitude  of  the  entire  signal  is  decreased  in  proportion  to  the  distance 
the  sound  traveled  from  the  source  to  the  receiver.  Second,  the  amplitude 
of  the  higher  frequencies  are  normally  diminished  more  than  the  lower  fre¬ 
quencies  because  of  the  attenuation  of  the  atmosphere.  The  addition  of 
background  noise  near  the  receiver  prevents  the  detection  of  weak  signals 
from  the  source.  Finally,  the  microphone  itself  often  introduces  some  dis¬ 
tortion  to  the  signal  due  to  a  phase  change  and  minimum  amplitude  cutoff 
caused  by  transforming  the  pressure  wave  in  the  air  into  an  electrical 
response  in  the  microphone.  This  imposes  some  analog  frequency  filter¬ 
ing  (often  to  avoid  aliasing  the  signal),  Hnally  transforming  the  analog  sig¬ 
nal  into  a  digital  signal. 


FFT  Mathematical  Output  of  Computer 

The  FFT  is  an  algorithm  that  decomposes  an  input  signal  reported  as  an 
amplitude  at  equal  time  increments  (referred  to  as  the  “time  domain”)  into 
both  an  amplitude  and  phase  at  equal  increments  of  frequency  (“frequency 
domain”).  The  transformation  assigns  the  largest  amplitudes  in  the  fre¬ 
quency  domain  to  those  frequencies  that  are  most  evident  in  the  time 
domain. 

The  discrete  FFT  is  deHned  as  follows:  given  amplitudes  g^  recorded 
at  time  t^  where 

tk  =  to  +  k  ■  ht  (2) 


and 

k  =  1,  2 . 2", 

then  the  FFT  defines  a  complex  number  F(fj)  at  each  real  frequency  fj 
fj  =  j  ■  if  ;  =  0,  1,  2,  ...,  2"-^  (3) 

where  the  frequency  increment  f  is  defined  as 
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(4) 


and  F(fj)  is  deflned  by  the  formula 

2-  (5) 

F{fj)  =  6r  •  5^  gjt  •  exp  (2kj  ■  ■  fj) 

h=l 

The  amplitude  A(fj)  and  phase  (fj)  are  defined  as  real  numbers  in  the 
normal  complex  variable  sense  (Ahlfors  1966)  so  that 

F(fp  =  A{fj)  exp  [2Ki  •  (6) 

The  amplitude  represents  the  strength  or  amount  of  that  frequency  fj 
present  in  the  time  domain  signal.  The  phase  represents  the  time  delay  or 
shift  from  a  cosine  wave  of  that  frequency. 

The  last  frequency  fj,  where  is  called  the  Nyquist  frequency.  It 

is  important  that  all  coherent  energy  for  frequencies  larger  than  the  Ny¬ 
quist  be  removed  before  computing  the  FFT.  Otherwise,  that  energy  will 
be  reported  incorrectly  by  the  FFT  as  some  energy  less  than  the  Nyquist. 
This  incorrect  reporting  is  called  “aliasing." 

All  of  the  amplitudes  and  phases  are  considered  to  apply  to  the  time  in¬ 
terval  over  which  the  summation  occurred:  [to,to+2“-6t].  This  frequency 
analysis  can  easily  be  repeated  for  later  time  intervals.  Figure  1 1  in  the 
main  text  shows  the  amplitudes  of  the  FFT  analysis  of  the  measured  data 
for  the  Hrst  sixty  frequencies  (at  a  1-Hz  increment)  and  for  one  hundred 
and  thirty-Hve  1 -sec-long  time  intervals.  That  graph  contains  three  bands 
of  Doppler-shifted  energy:  the  fundamental  of  the  main  rotor  blades 
(basic  frequency  at  17.2  Hz,  approaching  at  19  Hz  and  retreating  at 
16  Hz),  the  second  harmonic  of  the  main  rotor  blades  (basic  frequency  at 
34.4  Hz,  approaching  at  38  Hz  and  retreating  at  31  Hz),  and  the  third  har¬ 
monic  of  the  main  rotor  blades  (basic  frequency  at  SI. 6  Hz,  approaching 
at  57  Hz  and  retreating  at  47  Hz).  The  decrease  in  amplitude  at  the  higher 
harmonics  is  quite  evident. 

The  tabular  representations  of  parts  of  Figure  11  in  the  main  text  are 
shown  in  Tables  S,  8,  and  1 1  (again  in  the  main  text).  The  corresponding 
phase  information  appears  in  Tables  6,  9,  and  1 1  in  the  main  text.  The 
prominent  energy  of  the  fundamental  energy  appears  as  the  frequencies 
with  large  amplitudes  in  the  frequency  domain  that  can  be  tracked  from 
one  time  interval  to  the  next  (i.e.,  the  bold  numbers). 

The  main  restriction  of  the  FFT  is  the  restriction  that  the  discrete  steps 
5f  forces  on  the  frequency  resolution  of  the  energy  band.  The  main  pur¬ 
pose  of  the  Accurate  Calculation  of  Frequency  section  of  Chapter  2  of  this 
report  (pp  10-11)  is  to  obtain  an  improved  accuracy  over  the  normal  FFT 
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resolution.  This  higher  resolution  was  necessary  for  the  formula  for  the 
range  at  CPA  calculated  using  the  Doppler  theory  of  physics. 
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table  B1 

Description  of  Different  Frequency  Values _ 

Machwlcal _ 

■  Rate  of  revolution  of  moving  parts 
Aeouatlcally  Generated  at  Source  (1^.,  basic  frequencies) 

B  Restrict  mechanical  frequencies  to  those  that  are  audible 
B  Add  the  higher  harmonics  (see  Appendix  A) 

Acoustlcslly  Measured  at  Receiver 
B  Add  Doppler  shift  to  acoustically  generated  frequencies 
B  Deaease  amplitude  due  to  distance  traveled 
B  Impose  minimum  detectable  amplitude  due  to 
B  Limits  of  measuring  device 
B  Background  noise  exterior  to  measuring  device 
FFT  Mathematical  Output  of  Computer 

B  Report  amplitude  and  phase  at  equal  increments  of  frequency  for  any  signal  recorded 
at  equal  time  increments 

B  Impose  FFT  limitations 

B  Frequency  starts  at  0  Hz 

B  Frequency  increment  is  inverse  of  total  time  interval  analyzed 
B  Maximum  frequency  reported  is  Nyquist  frequency 

B  Frequencies  higher  than  the  Nyquist  that  are  present  in  the  input  time  signal  are 
irKorrecdy  reported  as  frequencies  below  the  Nyquist  in  the  frequency  domain 
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Table  B2 

Example  of  Mechanical  Frequency  Rates  for  Helicopters 


■  Main  rotor  shaft  revolution  rate  r^ 

{e.g.,  4.3  Hz) 

■  Main  rotor  blades  revolution  rate 

rg  m  (#  blades  on  main  rotor  shaft)  *  r^ 
(e.g.,  17.2  Hz  -  4  *  4.3  Hz) 

■  Tail  rotor  shaft  revolution  rate  ra 
(e.g.,  19.8  Hz) 

■  Tail  rotor  blades  revolution  rate 

r4  >  (#  blades  on  tail  rotor  shaft)  *  ra 
(e.g.,  79.2  Hz -4  *19.8  Hz) 

■  Turbine  engine  shaft  revolution  rate  rs 
(e.g.,  30C  Hz) 

■  Turbine  engine  blades  revolution  rate 

rg  >  (#  blades  on  turbine  engine  shaft)  *  rg 
(a.g.,  3,000  Hz  -  1 0  *  300  Hz) 


Table  B3 

Example  of  Acoustically  Generated  Frequency  Rates  for 
Helicopter  Sources 

Components: 

a  Main  rotor  blades  and  all  harmonics 
fih  ■  k  *  ra  for  k  ■  1 ,  2,  3, . .  . 

(e.g.,  17.2  Hz,  34.4  Hz,  51.6  Hz,  68.8  Hz.  86.0  Hz _ ) 

B  Tail  rotor  blades  revolution  rate 
fa  -  k  *  r4  for  k  ■  1,  2,  3, . .  . 

(e.g.,  79.2  Hz.  158.4  Hz.  227.6  Hz _ ) 

B  Turbine  engine  blades  revolution  rate 
fai,  •  k  *  rg  for  k  >1 1 ,  2,  3, . . . 

(e.g.,  3,000  Hz.  6,000  Hz.  9,000  Hz _ ) 

Composite  Results: 

Sequence  of  isolated  bands  of  energy  at  discrete  frequencies  (Hz):  1 7.2,  34.4,  51 .6, 
68.8,  79.2.  86.0.  103.2,  120.4,  137.6, 154.8,  158.4,  172.0.  189.2, .  .  . 
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Appendix  C 

Additional  Statistics  of  Other 
Signals 


All  real  data  possess  some  noise  that  causes  problems  with  well- 
intentioned  algorithms.  The  acoustic  data  of  Run  21  (Figure  11  in  the 
main  text)  is  no  exception.  Fortunately,  signals  other  than  the  fundamen¬ 
tal  of  the  main  rotor  can  be  used  to  estimate  some  of  the  information  con¬ 
cerning  the  helicopter.  Both  the  second  harmonic  (approaching  at 
approximately  38  Hz)  and  the  tail  rotor  (at  roughly  80  Hz)  can  be  tracked 
long  enough  to  make  some  of  these  measurements. 

These  are  essentially  independent  measurements  of  the  velocity,  range, 
and  the  time  of  the  closest  point  of  approach  (CPA).  The  variations  of  the 
separate  statistics  hint  at  the  magnitude  of  the  error  in  the  calculation. 
Since  the  input  data  to  the  calculations  are  real,  the  output  data  represent 
the  accuracy  of  both  the  algorithm  and  its  sensitivity  to  recording  and 
digitizing  noise.  Results  of  the  frequency  calculations  are  listed  below  in 
Table  Cl.  The  velocity  calculation  for  all  three  energy  packages  is  very 
similar  (see  Table  C2).  All  three  cross-spectral  density  (CSD)-derived 
velocities  are  above  the  velocity  reported  by  the  pilot  and  below  the 
velocity  predicted  by  the  Globtd  Positioning  System  (GPS).  The  CSD 


Table  Cl 

Summary  of  Frequency  Analysis,  All  Three  Strong  Signals 

Acoustic  Predictions 

Main  Rotor 

Tall  Rotor 

Main  Frequency 

2nd  Harmonic 

Approaching  frequency  (Hz) 

19.10 

38.20 

68.09 

Retreating  frequency  (Hz) 

15.65 

31.26 

72.20 

Basic  Frequency  (Hz) 

17.21 

34.38 

79.36 

Krwwn  frequencies  (Hz) 

17.2 

34.4 

79.3 

AppandIxC  AdcMionat  Statistics  of  Other  Signals 
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Table  C2 

Summary  of  Velocity,  Time,  &  Range  at  CPA,  All  Three  Strong 
Signals 

Acoustic 

Predictions 

Main  Rotor 

Tail  Rotor 

Pilots 

Reported 

GPS 

Estimates 

Main 

Frequency 

2nd  Harm 

Velocity  (m/sec) 

34.98 

35.20 

34.95 

33.46 

35.95’ 

Time  of  CPA  (sec) 

71.07 

71.39 

71.40 

61* 

3 

Range  at  CPA  (m) 

167 

166 

141 

182 

169 

Approach  velocity  average  of  time  intervals  (15,45]  A  [90,105]. 

‘  Indudas  1  .S'sec  addition  for  the  helicopter  moving  from  opposite  between  Arrays  A  and  B 
to  opposite  Array  A  and  0.5  sac  for  the  sound  to  travel  from  the  helicopter  to  the  receiver. 

^  There  was  no  independent  correlation  of  digitizing  time  to  GPS  time. 

velocities  must  therefore  be  considered  accurate  within  the  level  of  ac¬ 
curacy  for  the  independent  measurements  of  the  pilot’s  report  and  the  GPS. 

Ranges  predicted  at  various  times  for  all  three  signals  consistently  indi¬ 
cate  a  number  slightly  below  that  reported  by  the  pilots  (see  Figure  Cl). 
The  CSD  range  predictions  agree  well  with  Ac  estimate  of  the  GPS 
system. 
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Figure  C1 .  Best  three  signals’  prediction  of  range  at  CPA  utilizing  the  slope  of  Doppler  shift  -  Run  21  -  channel  1 .  Total  perpen¬ 
dicular  offset  =  1Q2  m.  The  pilot-reported  range  includes  the  0.52-sec  shift  from  CPA  to  receiver 
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A  Amplitude 

C  Velocity  of  sound,  m/sec 

CPA  Closest  point  of  approach 

CSD  Cross-spectral  density 

f(t)  Frequency  measured  at  receiver  at  time  t,  Hz 
f,  Far-field  approaching  frequency,  Hz 

fg  Basic  acoustic  frequency,  Hz 

ff  Far-field  retreating  frequency,  Hz 
FFT  Fast  Fourier  Transform 

g(t)  Fourier  approximation  of  h(t)  at  one  frequency 
GPS  Global  Positioning  System 

h(t)  Simulation  of  an  acoustic  amplitude  at  time  t 
//n(f)  FFT  over  [n,n+l]  of  the  signal  h(t)  at  f  hertz 
i  Complex  imaginary  unit  (-1)^ 

t  Time,  sec 

tcpy^  The  time  the  helicopter  was  at  CPA,  sec 
t(;pA(t)  Estimate  made  at  time  t  of  t^i^pA.  sec 

tg  Time  when  the  ground  is  receiving  sound  of  CPA,  sec 

AfipwNixD  Notation 


to(t) 

Estimate  made  at  time  t  of  t^,  sec 

R 

Range  at  CPA,  distance  of  helicopter  to  microphone,  m 

R(t) 

Estimate  made  at  time  t  of  the  range  at  CPA,  m 

TemCoh(f) 

Temporal  coherence  of  signal  at  the  frequency  f 

V 

Velocity  of  the  source,  m/sec 

WAM 

Wide-Area  Mine 

X 

Point  in  the  complex  number  plane 

^n(0 

TheCSD  of/fn(0 

e 

Phase,  degrees,  radians,  or  cycles 

P 

Phase  difference,  degrees,  radians,  or  cycles 

5f 

Frequency  increment,  Hz 

Izl 

Amplitude  of  the  complex  number  z 

* 

Complex  conjugate  operator 

D2 
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